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ABSTRACT 


This  paper  presents  the  results  of  an  analytical  and  experi¬ 
mental  study  of  thin  conical  and  hemispherical  shells  subjected  to 
combined  mechanical  and  thermal  loads. 

The  analytical  solutions  are  applicable  for  shells  constructed 
from  a  linear  viscoelastic  material  with  temperature  dependent 
properties. 

Experimental  models  were  subjected  to  various  constant  radiant 
heating  rates  to  various  steady- state  temperatures.  In  some  cases 
a  constant  normal  pressure  was  combined  with  the  thermal  load. 

Theoretical  values  for  the  meridional  and  circumferential  stress 
distributions  based  on  a  viscoelastic  analysis  (both  temperature  in¬ 
dependent  and  temperature  dependent  material  properties)  and  the 
elastic  analysis  are  compared  with  experimental  results  for  typical 
time  interval s  during  the  transient  and  steady-state  periods. 
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I.  INTRODUCTION 


In  recent  years  there  has  been  considerable  interest  in  the  class 
of  aircraft  problems  commonly  known  as  aerodynamic  heating.  The  effects 
of  aerodynamic  heating  on  the  structure  are:  (1)  thermal  stresses 
induced  in  the  structure,  (2)  creep  and/or  relaxation  effects,  and 
(3)  development  of  inhomogenity  since  the  material  properties  vary  with 
the  temperature  distribution. 

At  elevated  temperatures,  the  thermal  stress  problem  (Condition  1) 
is  but  the  first  of  several  related  problems.  The  phenomenon  of  creep 
(Condition  2)  becomes  of  importance  if  the  temperatures  are  high  enough, 
even  if  the  exposure  times  are  short.  In  addition,  creep  buckling  can 
take  place  since  the  stresses  are  compressive  as  well  as  tensile. 

to  * 

The  deformational  response  of  structural  metals  at  elevated 
temperatures  is  a  combination  of  elastic,  viscous,  and  plastic  compo¬ 
nents.  The  three  mechanical  constants  associated  with  these  compo¬ 
nents  of  deformation:  the  elastic  modulus,  the  coefficient  of  visco¬ 
sity,  and  the  yield  stress;  and  the  tkermal  constants  which  govern  the 
level  of  the  thermal  stresses:  the  coefficients  of  linear  expansion  and 
of  thermal  conductivity,  are  temperature  dependent  (Condition  3) . 

Classical  thermoelasticity  does  not  take  Conditions  2  and  3  into 
account  and  therefore  gives  only  crude  approximations  at  elevated  tem¬ 
peratures  when  these  effects  are  known  to  be  pronounced.  A  viscoelastic 
material*-  will  satisfy  Condition  2,  and  if  in  addition  its  material 
properties  are  taken  as  temperature  dependent,  then  true  aerodynamic 
heating  effects  are  approached.  A  viscoelastic  response,  that  is,  one 
which  exhibits  relaxation  and  heredity  effects,  can  be  assumed  to  repre¬ 
sent  approximately  the  behavior  of  metals  either  at  relatively,  low 
stresses  or  at  very  high  temperatures.  If  this  viscoelastic  response  is 
produced  by  temperatures,  the  term  thermoviscoelasticity  has  been  used. 

It  includes  the  general  theory  of  heat  conduction,  thermal  stresses,  and 
deformations  produced  by  thermal  flow  in  the  viscoelastic  media  and  the 
reverse  effect  of  temperature  distribution  produced  by  the  elastic  and 
viscous  deformation. 

The  first  concern,  therefore,  is  with  the  heat  conduction  in  the 
structure  itself.  The  usual  procedure  of  thermcelastic  analysis  is  to 
determine  the  stresses  on  the  assumption  that  the  thermal  and  the 
elastic  fields  are  uncoupled  and  that  the  temperature  field  is  given. 
Although  the  assumption  of  uncoupled  thermal  and  elastic  fields  has 
practical  truth  in  structural  problems,  the  two  fields  are  coupled  by  the 


Manuscript  was  released  by  the  authors  on  Nov.  I960,  for  publication 
as  a  ARF  Technical  Report. 

Numbers  refer  to  entries  in  the  bibliography.  The  bibliography  in¬ 
cludes  references  for  background  information  in  addition  to  those  references 
noted  in  the  text.  , 


second  law  of  thermodynamics.  The  heat  conduction  equation  will  contain 
elastic  terms  and  the  equations  of  elasticity  will  contain  temperature 
gradients.  The  conventional  formulation  in  thermoelasticity  simply 
assumes  that  the  thermal  field  does  not  contain  the  small  elastic  terms 
but  that  the  temperature  gradients  are  included  in  the  equations  of 
elasticity. 

The  equations  of  the  classical  theory  of  elasticity  may  be  immediat¬ 
ely  extended  to  viscoelasticity  "by  simply  replacing  the  elastic  moduli 
by  their  corresponding  operators?  As  in  the  case  of  thermoelasticity, 
the  reciprocal  coupling  between  the  temperature  and  deformations  is 
neglected.  Classical  thermodynamics  shows  that  one  effect  cannot  occur 
without  the  other  but  in  the  theory  of  structures,  its  order  of  magni¬ 
tude  is  not  significant. 

The  objective  of  the  program  reported  herein  which  has  been  con¬ 
ducted  at  the  University  of  Minnesota  under  Contract  No.  AF  33(616)- 
5723,  Task  70524,  Project  7063,  can  be  stated  in  general  terms  as 
"Development  of  Analytical  Methods  to  Predict  the  Structural  Behavior 
of  Curved  Surfaces  Subject  to  an  Environment  that  Simultaneously  Im¬ 
poses  High  Heating  Rates,  Temperatures,  Loading  and  Accelerations". 

Specifically,  the  investigation  is  planned  to  study  the  stress 
and  strain  distributions,  and  buckling  characteristics  of  surfaces  of 
revolution  subjected  to  external  pressure  and  various  temperature 
conditions. 

The  methods  of  analysis  take  into  account  the  influence  of  time 
and  temperature  upon  surfaces  of  revolution  with  temperature  dependent 
properties.  Conical  and  hemispherical  shells  are  selected  for  the 
verification  specimens  for  analytical  and  experimental  simplicity. 

The  methods  of  analysis  and  their  verification  by  experiments 
are  summarized  in  the  body  of  this  report.  Items  of  a  secondary  nature 
such  as  details  of  the  experimental  program,  are  presented  in  the 
appendices. 


II.  FUNDAMENTAL  EQUATIONS  OF  THERMOVISCOELASTICITY 


The  mathematical  formulation  of  the  response  of  a  viscoelastic  body 
to  combined  mechanical  and  thermal  loads  can  be  approximately  represent¬ 
ed  by  the  following  set  of  differential  equations: 

Equations  of  Equilibrium  »  o 

Compatibility  Equations  - 

Equations  of  State  P[ Sij  J 

2. 


where  P  and  Q  are  linear  operators  representing  different, viscoelastic 
materials  and  £ij  -  €ij -fa <fij  ,  - <JTj ~/3  ,  £ij  =&(1i,j+'fj,i) 

Subject  to  the  boundary  conditions,  the  system  of  nine  equations  on  the 
set  of  nine  unknown  field  variables  <Jcj }%i  is  complete  in  the  sense  that 
solution  of  the  system  is  unique,  if  the  solution  exists. 


2.1  The  Model  and  the  Equations  of  Equilibrium 

The  model  of  this  analysis,  as  shown  in  Fig  2.1-1,  is  a  thin  shell 
of  revolution  generated  by  any  arbitrary  curve.  By  thin,  it  is  implied 
that  the  shell  thickness  is  of  small  order  of  magnitude  as  compared 
with  other  dimensions.  The  shell  is  further  assumed  to  be  of  uniform 
thickness.  No  restrictions  are  placed  on  the  manner  by  which  the  edge 
is  supported.  However,  the  load  system  is  assumed  to  be  axial -symmetric, 
having  both  mechanical  and  thermal  origins. 

For  a  shell  of  revolution,  the  shear  stress  tangent  to  a  parallel 
circle  is  equal  to  zero,  or Qg  -  0s  .  Furthermore,  the  moment  and  the 
tensile  stress  in  9  -  direction  will  be  independent  of  9.  The  stresses 
and  moments  on  a  small  element  abed  defined  by  two  meridian  planes 
and  two  paralleled  circles,  are  shown  in  Fig  2.1  -  1. 

Three  equations  of  equilibrium  can  be  established  by  equating  the 
force  components  in  the  direction  tangent  to  the  meridian  and  normal  to 
meridian  to  zero  and  by  setting  the  moment  about  cd  to  zero. 

0  (2.1  -  1) 

+  (2.1  .  2) 

<P/V/  *  °  (2.1  -  3) 

Partial  differentiations  are  being  used  because  the  independent 
variables  are  also  functions  of  time. 

In  addition  to  the  three  equations  shown  above,  one  may  also  estab¬ 
lish  another  equilibrium  condition  by  considering  the  equilibrium  of  the 
portion  of  the  shell  above  a  parallel  circle 

2.THfysi"t+2nnQf>  coifl-hTVLp^o 

from  which 

/  ^  /  Pfl 

(2.i-4) 

One  way  to  simplify  the  equations  of  equilibrium  is  to  use  equations 
(2.1-2^,  (2.1  -  3),  and  (2.1  -  4)  and  intrprtucing  new  variables 
and  &s  described  in  Timoshenko2  ,  or  to  preceed  by  eliminating 

Qa*  from  the /equations  (2.1  -1),  (2.1  -  2),  and  (2.1  -  3)  to  obtain 
two  fundamental  equations  of  equilibrium.  This  procedure  of  approach 
had  been  used  by  Stodola23,  Keller *7 1  and  Frankhauser 2  in  their  elastic 
analysis.  A  brief  derivation  will  be  as  follows: 
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From  equation  (2.1  - 

Differentiation  with 


3)  one  may  write: 

Q  z  *  -M  cos  £ 

y  f,*t  4  r 

respect  to  4  yields 


(2.1  -  5) 

(2.1  -  6) 


Introducing  equations  (2.1  -  5)  and  (2.1  -  6)  into  equations  (2.1  -  1) 
and  (2.1  -  2)  gives 


(2.1  -  7) 


°{2.1  -  8) 


Equations  (2.1  -  7)  and  (2.1  -  8)  are  now  the  equations  of  equilibrium. 


2.2  Equations  of  State 


As  this  analysis  deals  with  linear  viscoelastic  material,  the 
relation  between  the  stress  and  strain  may  be  taken  as 

{S‘jj  (2.2-1) 


Upon  introduction  of  the  stress  and  strain  deviators, 
(2.2  -  1)  take  the  following  form 


=  *  (?£+<?)(*» -*T) 


equations 


(2.2  -  2) 


The  components  of  stress  may  be  written  as 

(jli+(i)(2ett  +  €«e~5e‘T)  1 

' =2\l7f+<*)[2£$9  +  fy-3tiT)  J  (2.2  -  3) 

Obviously,  both  <5^  and  CTe6  are  functions  of  the  physical  properties, 
strains,  time  anti  temperature.  It  will  be  discussed  later  in  paragraph 
2.6  that  the  physical  properties  involved  here  are  functions  of  temper¬ 
ature.  Since  the  temperature  may  vary  across. the  shell  thickness  as 
well  as  along  the  surface  meridian,  and  since,  the  deformations  for 
viscoelastic  material  are  functions  of  load  and  time,  the  resulting 
stress  pattern  is  extremely  complicated.  It  is  therefore  necessary  to 
make  reasonable  assumptions  to  simplify  the  problem  as  the  analysis 
proceeds. 


With  values  of  and  Gee  it  is  possible  to  write  the  expression 
for  membrane  forces  and  moments.  Taking  the  shell  middle  surface  as  a 
reference  surface  for  integration,  one  may  express  the  :..embrane  forcerf 
and  moments  as  follows 


A 


i{ 


,  Ma  f(°*>  ^('~k)Jr 


(2.2  -  4) 


But  for  thin  shells,  the  small  quantities  X  and  X  can  be  neg¬ 
lected.  Thus,  *f  A 


(2.2  -  5) 


These  integrations  and  the  values  of  ^  and  &qq  in  equation  (2.2  -  3) 
reveal  the  importance  of  simplification  in  all  variables  involved  in 
these  equations. 


2. 3  Strains  as  Functions  of  Displacement 


Considering  the  shell  middle  surface,  the  strain  or  deformation 
along  the  meridian  is 

(2-3-1) 

and  the  corresponding  strain  in  the  circumferential  direction  of  the 
parallel  circle  is 

^  =  A  ~M)  (2.3-2) 

Partial  derivatives  are  being  used  in  (2.3  -  1)  and  (2.3  -  2)  instead 
of  ordinary  derivatives  in  order  to  consider  displacements  as  function 
of  both  and  time.  If  the  variation  of  strains  across  the  shell 
thickness  has  to  be  considered,  th^strain  will  be  a  function  of  distance, 
f  ,  from  the  shell  middle  surface. 


(2.3  -  3) 
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Obviously,  the  strains  in  the  form  of  equations  (2.3  -  3)  will  add 
complication  in  the  derivation  of  basic  differential  equations.  In  the 
treatment  of  thin  shells,  however,  the  variation  of  deformation  across 
the  shell  is  neglected.  Appendix  AI  also  shows  the  justification  of  this 
approximation  by  the  order  of  magnitude  consideration.  With  this  approx¬ 
imation,  equation  (2.2  -  3)  may  be  written  . 

(2  ^6/1  +  J  (2.3-4) 

In  other  words,  in  these  stress  equations,  we  take  the  strain  of  the 
shell  middle -surf ace  as  the  mean  value  of  the  strain  across  the  shell 
thickness. 


2.4  Pressure  and  Temperature  of  the  Shell 

At  any  point,  the  shall  is  subjected  to  a  pressure  load  p,  which 
could  be  taken  as  a  function  of  time  or  a  constant  depending  upon  the 
status  of  the  loading,  steady  or  transit.  In  the  case  of  a  flight 
vehicle,  the  pressure  is  of  dynamic  type.  The  boundary  layer  immediately 

outside  the  body  surface  and  the  phenomena  of  separation6  result  in  a 
variation  of  pressure  along  the  meridian  of  a  shell.  The  pressure  is 
maximum  at  the  stagnation  point  and  varies  as  one  moves  downstream  on 
the  body.  Furthermore,  if  the  angle  of  attack  is  not  zero,  the  pressure 
distribution  will  be  unsymmetrical.  However,  this  analysis  is  restricted 
to  bodies  of  zero  angle  of  attack  so  that  the  pressure  load  will  be 
axial-symmetric.  Fig  2.4  -  1  shows  the  pressure  distribution  along  a 
meridian. 

In  the  case  of  high  speed  vehicles,  aerodynamic  heating  effects  are 
experienced  in  addition  to  the  surface  pressure  distribution.  The 
resulting  surface  temperature  distribution  normally  varies  along  the 
body  surface.  The  manner  of  variation  depends  upon  the  flight  Mach 
number  and  the  heat  transfer  conditions.  In  this  analysis,  the  body  has 
a  constant  inner  surface  temperature,  with  heat  conduction  on  the  edge. 
Thus  the  resulting  outer  surface  temperature  will  be  shown  as  in  Fig 
2.4  -  2.  In  some  cases,  the  mathematical  complication  arising  from  the 
variation  of  temperature  along  the  meridian  forces  one  to  assume  a  uniform 
outer  surface  temperature.  This  assumption  deviates  quite  seriously 
from  the  actual  problem.  However,  excellent  approximation  can  be 
made  by  taking  a  number  of  uniform  temperatures  as  shown  in  Fig  2.4  -  3. 

This  analysis  is  primarily  concerned  with  the  temperature  variation 
within  the  shell  thickness.  For  steady  state  condition,  the  temperature 
gradient  across  the  shell  thickness  is  a  constant  (based  on  the  assump¬ 
tion  of  constant  inner  surface  temperature  and  thickness).  For  the 
transient  state,  the  temperature  gradient  varies  with  the  outer  surface 
temperature  and  with  the  distance, from  the  middle- surface. 
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2.5  Boundary  Conditions 

With  the  pressure  and  temperature  prescribed  in  Paragraph  2.4, 
the  boundary  conditions  can  be  derived.  Depending  on  the  choice  of  the 
initial  conditions,  one  of  the  following  catagories  of  boundary  con¬ 
ditions  should  apply: 

a)  Stress  boundary  conditions, 

b)  Strain  boundary  conditions, 

c)  Temperature  boundary  conditions, 

d)  Displacement  and  rotational  boundary  conditions. 

It  is  not  intended  that  a  detailed  discussion  of  all  possible 
boundary  conditions  be  presented  here.  However,  in  the  subsequent  me¬ 
thods  of  analysis,  selection  of  appropriate  boundary  conditions  are 
illustrated. 

For  the  problem  under  discussion,  the  following  diagram  summarizes 
the  principal  boundary  conditions  applicable  for  desired  solutions: 


(see  page  8) 
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BOUNDARY  CONDITIONS  DIAGRAM 


Fixed- supported  Edge 


Simple -supported  Edge 


At  t  =  0,  the  shell  is 
at  uniform  temperature. 
The  initial  stress  and 
strain  will  he  that  re¬ 
sulting  from  the  pres¬ 
sure  load  alone 


At  t  =  0,  a  steady  tem¬ 
perature  gradient  exists 
across  the  shell  thick¬ 
ness.  The  initial  stress 
and  strain  will  be  that 
resulting  from  both  pres¬ 
sure  and  thermal  loads 


At  t=0,  the  tempera¬ 
ture  gradient  across 
the  shell  thickness  is 
a  function  of  time. 

The  initial  stress  and 
strain  will  be  that 
resulting  from  both 
pressure  and  thermal 
loads  corresponding 
to  t«=0 


At  t  =  0 


At  t 


0 


p  =  0 


p  4  o 


At  t  =  0 


At  t  =  0 


Material  is  Elastic  Material  is  Viscoelastic 


2.6  Physical  Properties 

The  physical  properties  of  shell  material  involved  in  this  analysis 
are:  viscosity,  shear  modulus  of  elasticity  and  thermal  expansion  co¬ 
efficient.  For  the  most  metals,  viscosity  and  shear  modulus  of  elas¬ 
ticity  decrease  with  an  increase  in  temperature,  while  the  thermal  ex¬ 
pansion  coefficient  increases  with  temperature.  Since  we  are  concerned 
with  the  effects  of  aerodynamic  heating  on  viscoelastic  bodies,  the 
variation  of  material  properties  with  temperature  has  to  be  taken  into 
consideration. 

Fig  2.6-1  (Reference  4)  shows  the  variation  of  the  mean  coefficient 
of  thermal  expansion  with  temperature  for  several  typical  materials.  It 


"V 
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s. 


is  obvious  a  linear  relationship  will  closely  approximate  the  actual 
variation.  However,  for  some  materials,  the  temperature  effect  on  the 
thermal  expansion  coefficient  is  small  as  compared  with  other  properties. 
For  example,  in  case  of  aluminum,  the  value  of  thermal  expansion  coef¬ 
ficient  is  increased  by  8%  for  a  temperature  variation  from  250°F  to 
600°F.  However,  for  tne  same  temperature  variation,  the  shear  modulus 
may  decrease  as  much  as  50%.  In  the  subsequent  analysis,  therefore,  the 
thermal  expansion  coefficient  is  assumed  to  be  constant. 

The  shear  modulus  of  elasticity  takes  a  similar  trend  of  variation 
with  temperature  as  the  modulus  of  elasticity.  The  value  of  shear 
modulus  of  elasticity  sometimes  can  be  approximated  by  faking  one  third 
that  of  modulus  of  elasticity.  The  variation  of  modulus  of  elasticity 
is  shown  in  Fig  2.6  -  24.  Similar  to  the  coefficient  of  thermal  expan¬ 
sion,  a  linear  variation  may  be  assumed.  The  following  relationship  is 
used 

=  $.(/- 

is  the  shear  modulus  at  some  datum  temperature 

is  the  shear  modulus  at  T  degrees  above  the  datum  tempera¬ 
ture 

is  the  average  rate  of  change  of  G  with  respect  to  the 
temperature  for  a  specific  temperature  range. 


where 

fi 


The  information  of  the  viscosity  for  metals  is  very  limited.  Most 
liquids  follow  the  Newtonian  Law,  that  is, the  viscosity  decreases  ex¬ 
ponentially  with  the  inverse  of  the  absolute  temperature.  One  may  use 
a  similar  function  for  metals  but  the  integration  of  viscosity  with  respect 
to  temperature  T  will  be  extremely  difficult.  To  retain  the  exponential 
property  of  viscosity  on  one  hand  and  to  take  the  advantage  of  less 
definite  viscosity  theory  of  metal  on  the  other  hand,  it  is  assumed  that 
the  viscosity  varies  according  to 

is  viscosity  at  some  datum  temperature 

is  viscosity  at  a  temperature  T  degrees  above  the  datum 
temperature 


where 


1° 

1 


/<  is  evaluated  between  the  datum  temperature  and  a  prespeci¬ 
fied  temperature. 


2.7  Methods  of  Analysis 

In  paragraph  2.1  it  is  shown  that  there  are  two  ways  of  combining 
the  equations  of  equilibrium  which  suggest  there  are  different  ways  of 
solving  the  problem.  Paragraph  2.3  indicates  that  the  strain-displace- 
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ment  relation  may  take  a  more  simplified  form.  The  physical  properties 
may  betaken  as  functions  of  temperature,  or  they  may  be  evaluated  at  some 
appropriate  mean  temperature  in  which  case  the  problem  is  treated  as 
being  independent  of  temperature.  All  of  these  choices  suggest  that 
there  are  many  ways  to  arrive  at  a  solution.  In  Section  III  the 
conical  shell  is  analyzed  by  three  different  approaches.  However, 
from  a  mathematical  point  of  view,  for  a  given  shell  geometry  one 
method  may  be  more  efficient  than  the  other. 

The  general  procedure  of  analysis  is: 

1.  To  procede  with  the  given  temperature  profile  which  leads  to 
the  physical  properties  and  thermal  stresses. 

2.  To  make  a  choice  of  strain-displacement  relation  for  substitu¬ 
tion  into  the  equation  of  state. 

3.  To  find  the  values  of  membrane  forces  and  moments  from  the 
stress  components  and  express  them  in  terms  of  displacements,  material 
properties,  and  independent  variables.  In  some  cases,  the  use  of 
appropriately  chosen  new  variables  may  be  useful. 

4.  To  introduce  the  values  of  membrane  forces  and  moments  into 
the  equations  of  equilibrium  to  obtain  the  governing  equations.  In  some 
cases,  instead  of  using  membrane  forces  as  obtained  from  3  a  fourth 
equation  of  equilibrium  may  be  used. 

5.  To  solve  the  governing  equations  the  boundary  conditions  are 
inserted  to  evaluate  arbitrary  constants  that  arise  from  the  integration 
of  the  differential  equations. 

The  process  of  analysis  may  be  summarized  with  the  following  diagram 


PROCESS  OF  ANALYSIS  DIAGRAM 


Temperature  Profile - 1 

Physical  Properties 


-Strain-Displacement  Relations 


EquationSjof  State 
Membrane  Forces  and  Moments 
Equations  of  Equilibrium 

I 

Governing  Eqi  „_ions 


•Boundary  Conditions 


Solutions 
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III.  ANALYSIS  OF  CONICAL  SHELL 


The  critical  section  of  many  space  flight  vehicles  is  of  conical 
geometry.  The  conical  shell  therefore,  is  taken  as  one  of  the  examples 
for  demonstrating  the  analysis. 

Three  different  methods  of  analysis  are  used  to  attack  this  problem 
and  solutions  are  illustrated  for  specific  conditions. 


3.1  Simplified  Strain-Displacement  Relation  Method 

3.1.1  The  Governing  Equations 

For  the  conical  shell,  it  is  convenient  to  use  the  shell  element 
length  1  as  an  independent  variable  which  is  measured  from  the  tip  of 
the  shell.  The  principal  radius  J>  approaches  infinity  with  approaching 
eU  and  fz=jrAuY,  Y  being  half  of  the  apex  angle. 

The  equations  of  equilibrium  are  now  written  as 


-?M<L  „  ye 
2/ 


=  o 


“>rY  +  o 


(3.1.1  -  1) 


9J 

which  may  be  combined  into  two  equations 

-  f/a 


_  21^1-+.^/  ccrY +  e{-=  O 

?r  m 


(3.1.1  -  la) 


For  the  strain-displacement  relations,  the  simplified  functions  dis¬ 
cussed  in  paragraph  2.3  and  also  shown  in  Appendix  Aiare  used 

*  -JT 

,  .  i/~-  us  car  y 

"  y 

A  uniform  outer  surface  temperature  is  assumed  to  be7t  degrees  above  the 
inner  surface  temperature.  However,  ^  may  be  considered  as  a  function 
of  time,  if  desired.  The  steady-state  temperature  gradient  across  the 
shell  thickness  will  be  {5 //) and  the  temperature  at  a  distance  %  from  the 
middle  surface  will  be  T  degrees  above  that  of  the  inner  surface  tem¬ 
perature,  or 


(3.1.1  -  2) 


The  temperature  dependent  properties  are  written  as 
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T-ie-^ 


(3.1.1  -  3) 


oi  =  <*„ 


Introducing  (3.1.1  -  3)  into  equation  (2.3  -  4),  one  finds  for  the 
stress  components 

f  !/  rfl  D  f~  /  /  f 

°i>r2 

*  2 


(3.1.1  -  4) 


Thus  the  membrane  forces  and  the  moments  can  be  integrated  as  required 
by  equations  (2.2  -5).  For  example 

-  -m 


r 

2  6 


— ^  ft  p 

-  *P,t  r" ¥+ik(z »f 


nt 


yj 


(3.1.1  -  5) 


M,  -  2/. « 

^  ' 


tilk-2 


■]  2H+- 


{Tt  j  3t 


-h 


G^TtS1 

2 


(3.1.1  -  6) 


Introducing  equations  (3.1.1  -  5)  and  (3.1.1  -  6)  and  the  corresponding 
expressions  for  Ma  and  into  the  equation  of  equilibrium  (3.1.1  -  la) 
one  obtains  r  r 
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(3.1.1  -  7) 


and 


21^  2  + 

+  torts,** &  aW»5*J?  blklLJp 

^  2  ^  /  ‘L  7. ~ 

-f  Z<£,J”  ^/-  or  v(2&6l1+  t  pf  - 

-  or^y.-  e'^4/**  ¥■  - 


_  lu»t4££- 


?)]*{*-?*) 


W£i 


(3.1.1  -  8) 


3,1.2  Method  of  Solving  the  Governing  Equations 


First  of  all,  introducing  new  variables  T T~2€fai-6an  and to 
replace  terms  involving  and  (=9/ut  ,  results  in  the  following  form 

for  the  governing  equations 


m)M*  *?+37Zv)-h  uh° 

and 

%t&) ^H-k 

+  &£gSLM^+,^+<&(2-fii)<.rr(iir-3r)+r4- 


(3.1.2  -  1) 


-c1°T^le  +-  Jz§~)*r  r-jJf-+6J*rt(3-2/iK)o’T r  (3.lt2  -  2 > 

From  equation  (3.1,2  -  1) 

and  for  smoothly  continuous  functions  of U  and  \T with  respect  to  /  and 
■£  ,  it  can  be  written  that 

ft  jf*SKlj 

Hence 
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,  fir 

'M£l 


+-  3 


2>lr 

?tjt 


Using  equations  (3.1.2  -  1), 
governing  equations  become 


(3.1.2 


(3.1.2  -  2),  and  (3.1.2  -  3),  the  new 


3) 


(3.1.2  -  la) 
(3.1.2  -  4) 


where 


r2  = 


f,- 


Z- 


7X.<s?T 

2 fy(i-e 

*Tt 


i+e 


!-c 


MZ. 


,-^r 

j 


a  function  of  time,  then  Fi ,  F2,  Fo,  and  Fa 

mi _ _ a_i _ a  _  n  _ a.  • _ _  n  _ •r3a_i _ _ *7 


VO-*-*! 

For  cases  where  Tt  is 

are  all  functions  of  time.  The  method  of  separation~of  variables  are 
used  to  solve  the  governing  equations.  By  trying  a  solution  ((3%i+3V-Tf\ 
of  equation  (3.1.2  -  la)  in  the  form  Ljf/Ti (ir)  ,  one  obtains 

r  ^ 

(e%<-sr-vj.  L,(tme  (3.1.2  -  5) 


The  values  of  Uf)Z(o)  may  be  obtained  from  the  boundary  conditions, 
i.e.,  the  value  of  ( £  -  V+-3V)^t  ^  at  a  point  Jl  and  at  an  ar¬ 

bitrarily  chosen  initial  time.'  It  is,  therefore,  the  choice  of  boundary 
conditions  which  has  an  influence  on  whether  the  solution  of  the 
problem  can  be  obtained.  For  this  example,  the  boundary  conditions  are 
chosen  such  that 


a)  A  uniform  and  constant  pressure  is  exerted  on  the  shell  at  all 
times. 

b)  It  is  a  function  of  time,  equal  to  zero,  at  i  =  0;  increasing 
to  a  constant  value  when  a  steady  state  condition  is  reached^  For 
simplification  it  is_assumed  that  7t  =  C-b  for  0 ^£^4  and  Tt-Tt  for 
t>tc  3  where  C  and  are  constants  and  £  is  the  transition  period. 

c)  The  edge  of  the  shell  is  fixed. 
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These  boundary  conditions  imply  that  at  =  0  the  shell  has  a  defor¬ 
mation  which  is  caused  by  the  pressure  load  alone.  The  next  question 
is  the  state  of  the  material  conditions  at^«*  0,  whether  elastic  or  vis¬ 
coelastic.  This  depends  on  the  type  of  material  and  the  temperature  to 
which  the  material  is  subjected  at  the  initial  time.  In  the  elastic 
case,  one  may  evaluate  (/ from  results  of  elastic  analysis 
as  discussed  in  most  applied  elasticity  textbooks.  A  simplified  version 
of  the  elastic  analysis  of  a  conical  shell  with  thermal  stress  is  dis¬ 
cussed  in  Appendix A2  .  For  the  present  application,  we  will  use  the 
results  of  AppendixA2  but  omit  the  contributions  from  thermal  effects. 
However,  if  the  boundary  conditions  are  specified  as  an  elastic  material 
with  thermal  stress  at  the  initial  moment,  then  the  complete  results 
should  be  used. 


For  elastic  analysis 


A., = £  fa*  K)a  ,  -  fs 


twY  1 


(3.1.2  -  6) 


A^  and  A2  are  Lo  be  evaluated  from  the  edge  conditions. 

The  equation  (3.1.2  -  5)  may  be  written  ^ ^ 

(3.1 . 2  -  7) 

From  (3.1.2  -  7),  the  values  of  ^%t  can  be  expressed  in  the  derivatives 
of  TT  and  the  known  initial  values  , 


if- 


where 


(3.1.2  -  8) 

±t>wY 


substitution  of  ^/l from  equation  (3.1.2  -  8)  into  equation  (3.1,2  -  4) 
reduces  it  to  a  single  dependent  variable  equation. 


C;d  1  L  J  EL 


(3.1.2  -  9) 
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Again  using  the  method  of  separation  of  variables  and  boundary  conditions, 
one  finds  for  ~U~ 

TT=  (  [ I \fc-Fj 


(3.1.2  -  10) 

From  equation  (3.1.2  -  7)  the  corresponding  value  of  V  is  found  to  be 

Jf&f  riff  frjt  ,  T 

r-  -£j-  (i[l  1 

+  0-2*)  2_(k;  A/}  \ 

Eel  91  W+'fyeJ 


(3.1.2  -  11) 


With  values  of  TT  and  V~  the  subsequent  values  of  '2-^n  +  £&r/  and 
2  F0n  +  FJh  their  time  derivatives,  and  the  membrane  forces  and 

moments  rcan  be  readily  determined.  The  values  of  fi/0  and  Ffg  are  illus¬ 
trated  as  follows  t 

-f FtJt  r 

'  -  ^ 


zfyO-e**),  F 

v-r  + 


Kii 


(3.1.2  -  12) 


and 
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(3.1.2  -  13) 


From  this  example,  it  can  be  seen  that  the  method  of  analysis  can 
be  extended  to  apply  to  other  boundary  conditions;  for  example,  if 
pressure  is  being  a  function  of  time  and  distance  along  the  shell  element, 
pl-t/i)  ,  or  if  an  initial  temperature  gradient  across  the  shell  thickness 
exists  or  if  the  edge  is  simply  supported. 


3.2  Simplified  Differential  Equation  Method 

One  of  the  important  simplification  made  in  paragraph  3.1.1  was  the 
omission  of  terms  containing  ^  in  the  equations  of  and  ^  .  In 

doing  so  it  was  assumed  that  the  middle  surface  values  were  the  mean 
values  and  that  the  order  of  magnitude  of  these  terms  was  small.  In  this 
paragraph,  the  basic  governing  equations  are  derived  without  imposing 
this  limitation,  however,  J  and  w  are  assumed  not  to  be  functions  of  ?  . 


write 


Applying  equation  (2.3  -  3)  to  the  conical  shell  results,  one  may 


r  ar’ 

_  if— uj  cor  IT  T  -3a/ 
1  /  zi 


(3.2  -  1) 


Employing  the  same  temperature  profile  as  given  in  paragraph  3.1.1,  the 
stress  components  may  be  written  as 


n 


-far)' 


if-  ajcorf  f  ^ko 

t  l  H 
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(3.2  -  2) 


<£-: ffc]l[*(eTcfid*#-t&  -*44) 


The  membrane  forces  and  moments  will  be 
.1 


-  us*ee—  [- 


A 

x-L^r 


**)*(* 


*2V('-*Kfr*  **?*)-  ^  &)- 

£  T/. 


(3.2  -  3) 


(3.2  -  4) 


(3.2  -  5) 
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(3.2  -  6) 


(3.2  -  7) 


Introducing  equations  (3.2  -  4)  through  (3.2  -  7)  into  the  equations  of 
equilibrium  (3.1.1  -  la),  one  reduces  the  governing  equations  to  two 
5th  order  linear  partial  differential  equations.  They  are 
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(3.2  -  8) 


+¥{44i(&*£-6£+££)+ 


**  r 


*  n  ^-kA «  f  2^}mj(,-^2^)^ 


(3.2  -  9) 

The  order  of  magnitude  of  each  terms  in  these  two  equations  is 
investigated  by  a  method  similar  to  that  discussed  in  Appendix  A.  With 
j  and  tv  having  the  order  of  magnitude  of  shell  thickness  /  ,  the  terms 
involving  ir  and  tv  in  equation  (3.2  -  8)  have  the  order  of  magnitude  £/%) 
'/M  ,  ft  ,  ™d  r/‘)‘  V  “d  i«  equation  (3.2  -  ?,(%)’  h 

s('/c) ,  t(*/ty  and//%;  • 

,  Neglecting  the  terms  having  the  order  of  magnitudes  oi'/t(/t)  and 
r/l)1  in  equation  (3.2  -  8)  and  %(*/)* ,/lyt)1  terms  in  equation  (3.2  -  9), 
the  governing  equations  are  written  as  follows 


*24.(1-  -  **?*)+  */&- »T r &j=. 


(3.2  -  10) 
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,  tt%. 

4  2<gj(,~  £Z)(  jf+2-~-*J*rr+-  k *££.  -Sarrfc&^^Jg^g- 

-<*><?<*%  corf  (j-2fTtJ+ [4=0 


(3.2  -  11) 


Comparing  equations  (3.2  -  10)  and  (3.2  -  11)  with  the  corresponding 
equations  (3.1.1  -  7)  and  (3.1.1  -  8)  which  were  derived  by  applying  the 
order  of  magnitude  consideration  before  the  derivation  of  the  governing 
equations,  one  finds  the  first  equations  are  exactly  the  same  while  the 
second  equations  differed  by  the  terns 


From  a  mathematical  point  of  view,  the  solutions  of  equations 
(3.2  -  10)  and  (3.2  -  11)  will  a  better  approximation  than  that  of 
equations  (3.1.1  -  7)  and  (3.1.1  -  8),  provided  that  both  sets  of  equa*= 
tions  can  be  solved  analytically.  However,  it  has  been  found  that  in 
this  case  the  governing  equations  are  not  as  readily  solvable  as  in  the 
former  case.  Some  efforts  have  been  spent  to  simplify  the  governing 
equations.  The  results  are  presented  here. 


(3.2  -  12) 


Zr~3Tr] 
'  / 


-A 


(3.2  -  13) 


where 


TT=  z-jf+-jr(w-  aJ“>rY') 
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£ 


=  2  / 


/+C 


■KK 


he 


-  *-n 


&o<f /Ik /t  *7)n/ ir 


and  f},  F,  ,  Ft 


are  defined  in  paragraph  3.1.2. 


3.3  Constant  Property  Method 

Although  the  main  interest  of  this  analysis  is  concerned  with  shells 
which  have  tempera£ure  dependent  material  properties,  a  crude  approxima¬ 
tion  however,  can  be  made  by  assigning  properties  at  some  arbitrary 
constant  temperature,  say  the  middle  surface  temperature,  and  then 
treat  the  problem  as  a  temperature-independent  one. 


3.3.1  Governing  Equations 

The  equations  of  equilibrium  will  be  the  same  as  those  in  paragraph 
3.1.1. 


2feklL=M, 

21 

cor  if  +  pf = o 


(3.3.1  -  1) 


(3.3.1  -  2) 


The  strain  components  remain  as 

<5r 

-  2 (l 2 

except,  here  1 ,6j/  and  oi.  are  constants.  Temperature  profile  across  the 
shell  thickness  is  assumed  to  be  the  same  as  given  in  paragraph  3.1.1, 
i.e.,  f~  (fa- ■ 


Using  the  strain-displacement  relations  as 

/l  _  u~-  uj  /of  S'  ?■ 

- t - T2t 

In  this  method,  the  procedure  may  follow  that  given  by  Timoshenko 


(3.3.1  -  3) 
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Intro due ing 


the  moments  can  be  written 
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or 


The  value  of  /i^  is  derived  from  equation  (2.1  -  4)  as  follows 

lEl&nir  fycor  IT-*-  2zltyss«*r+  rr(f#S  *)r° 


ty~~  fy’*"*'-  ^ZMT”~(^*'^)rMr  (3.3.1 

From  the  third  equilibrium  equation  of  (3.3.1  -  1),  Mg  is: 

(3.3.1 


-  5) 


6) 


Combination  of  equations  (3.3.1  -  4)  through  (3.3.1  -  6)  and  the  intro¬ 
duction  of  the  new  variable  W-ty! 7*u)f  yields  the  first  governing  differ¬ 
ential  equation  r 


-+&2L, 

31'  J  31' 


!  &  2~C  -L^sL.  ,.<Y—  3cotY W~ 3  et  T/t 

1  tag  ft  h  ~r  to  yr  *  z  *  m/ 


3U6  IT 
at 


(3.3.1  -  7) 


The  integrated  equation  of  state  in  paragraph  2.2 

A •  t* 


4' 


-  g  1 


7  ’ 


n-t 


may  be  used  to  derive  the  strain  components.  They  are 


and 


_  c 


-f* 


&}  J° 

Q,  -t 

-V  Ci 


r  cCt  +*/ 


4.  .  C_L  ++ / 


(3.3.1  -  8) 

(3.3.1  -  9) 

(3.3.1  -  10) 


The  new  variable  ^  may  be  also  written  in  terms  of  the  strain 
components. 


(3.3.1  -  11) 
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Combining  equations  (3.3.1  -  5),  (3.3.1  -  6),  and  (3.3.1  -  9)  through 
(3.3.1  -10),  one  may  write  the  second  governing  differential  equation 
as 


§Nr-f  - 


(3.3.1  -  12) 


Equations  (3.3.1  -  7)  and  (3.3.1  -  12)  are  now  the  governing  equations  of 
the  thermal  viscoelastic  conical  shell  of  temperature  independent  proper¬ 
ties. 


3.3.2  Method  of  Solving  and  Solution 

,  From  equation  (3.3.1  -  7)  the  variable  ~W may  be  expressed  in  terms 
of  sp  and  its  derivatives. 

W~~  f  Mi  1  H  (3.3.2  -  1) 

After  computing  derivatives  of  ~W~ with  respect  to  /  and  substituting  the 
results  into  equation  (3.3.1  -  12),  one  obtains  the  fifth  order  partial 
differential  equation 


(3.3.2  -  2) 


It  is  easily  seen  that  7t  in  this  method  of  analysis  has  much  more 
flexibility;  it  may  vary  with  the  time  t  as  well  as  with  the  shell 
element  length  /;  i.e.,  Ji  ~  7t  (4 ,i )  .  For  uniform  outer-surface  tem¬ 
perature  ^  ;  and  for  steady  state  thermal  load  ^7* /&(■  -  o  .  Method  of 

the  separation  of  variables  may  be  used  to  separate  equation  (3.3.2  -  2) 
into  two  ordinary  differential  equations;  a  first  order  equation  with  re¬ 
spect  to  time,  and  a  fourth  order  equation  with  respect  to  the  shell 
element  length.  The  complementary  equation  of  the  fourth  order  equation 
may  be  further  separated  into  two  Bessel  equations.  Their  solutions  are 
the  Bessel  functions  of  the  first  kind  and  the  second  kind  in  complex 
variables,  which  can  be  expressed  in  ber,  bei,  kei,  ker  functions.  With 
the  choice  7**'*'/'/]  as  a  particular  solution,  the  complete  solution 

for  uniform  and  steady  load  condition  may  be  readily  written  as  follows 
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I 


'f'  ^j*)+ -^-)*£r ”^[ie  ^ 26s\ 

K-2 y^^r[A(^-^fy^-hA(k  M~)]-r/a“r 
q- s*£«rr[Aff^)  “ 

Mt '  ~z£arri '[A(^+^~^f)+A(  )+t{-)+M"-)]-  ^2^ 


J 

(3.3.2 


where 


Avji 

?/?)*“  a***? 


J?'=#L 

7: F 


and 


z  /2  CorY~ 


3) 


7 

A  ,  A.,  A2,  A-,  and  A,  are  to  be  determined  by  the  boundary  conditions. 
A°  will  be  determinedly  the  load  condition  at  the  initial  time,  i.e., 
t°-  0.  The  continuity  of  ^at  ,/=<?  inmediately  eliminates  A^  and  A^. 

The  edge  conditions,  then,  determine  the  values  of  A^  and  A2. 


IV.  ANALYSIS  OF  HEMISPHERICAL  SHELL 


The  analysis  discussed  in  this  chapter  has  greater  applicability 
than  the  title  suggests.  It  may  be  applied  to  any  spherical  shell  having 
a  maximum  £>  equals  to  or  less  than  ufz  • 

The  method  used  here  is  similar  to  that  described  in  paragraph  3.1, 
i.e.,  the  simplified  strain-displacement  relation  method.  However,  it  is 
not  suggested  that  this  is  the  only  method,  on  the  contrary,  the  other 
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two  jnethods  presented  in  Paragraphs  3.2  and  3.3  may  be  applied.  However, 
the  application  will  not  be  discussed  to  avoid  repetition. 


4.1  Governing  Equations 

The  independent  variable  ^  used  in  the  equations  of  equilibrium 
(2.1  -  1  to  -  3)  could  be  retained.  The  values  of  J}  and  are  the  same 
and  are  equal  to  the  radius  of  the  sphere  fa  .  The  radius  of  parallel 
circles  ft  may  be  replaced  by  Pstud  .  Hence  the  equations  of  equil¬ 
ibrium  Jo  r 


^ -/^g  as?- o  (4.1  -i) 

fjfsvf-o  (4.1  -  2) 

(4.x.  3) 

To  el iminate  Qf  f„  (4.1  -  1)  is  differentiated  with  respect  to 

and  is  added  to  equation  (4.1  -  2).  Thus 


Subtraction  of  equation  (4.1  -  1)  from  equation  (4.1  -  3)  yields 


'"W  V  / 


7 S  +/>  fa-fy)  cor /=  * 


(4.1 


4) 


5) 


The  two  equations  of  equilibrium  (4.1  -  4  and  4.1  -  5)  are  then 
used.  Consequently,  the  expressions  for  f  AA#,  AAs  and  A/tf  are  the 
same  as  those  of  equations  (3.1.1  -  5)  through  (3.1.1  -  8).  Substitution 
of  the  latter  equations  into  equations  (4.1  -  4)  and  (4.1  -  5)  yields 
two  governing  differential  equations 


-I 2  ^ 

kJt  l  'L  3  r  ty*  Mr  j 

+  31)  jk&0al  *  3  J-z(z6>„+ - 

-  ***(?-*&* 

and 
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(4.1  -  7) 


where 


4.*zfps(fsB-2f.l 


F,~ 


(4.1  -  7a) 


4.2  Method  of  Solution 


Introduction  of  new  variables  ~TT-  2  fyh  *  and  I/=  (=pM  into  the 
governing  equations  (4.1  -  6)  and  (4.1  -  7)  yields 


and 


-  ~r) 

+  ?  /-fr  +  "r/frr-v] 


(4.2  -  1) 


(4.2  -  2) 


Using  the  method  of  the  separation  of  variables  in  equation  (4.2  -  2) 
one  finds  _  -fpj, 

2E.+  cr?{3r-o)=/j^+c»rtf3r-vjJ e  «' 

y  ^  ^  (4.2  -  3) 

The  value  of  (<tS  (3vr-°) Jt  .0has  to  be  evaluated  from  the  boundary 

conditions  at  the  initial  time,  i.e.,  the  temperature  profile,  the 
pressure  load,  the  properties  of  material  and  the  stress  and  strain 
induced  at  t  *  0.  If  it  is  an  elastic  shell  under  pressure  load  with¬ 
out  thermal  stress,  the  analytic  solution  of  elastic  hemispherical  shell 
given  in  various  mechanics  textbooks  can  be  used  to  evaluate  the  initial 
condition.  Thus  knowing 


/QiMnc 


and 


feJL- 
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By  the  definition  of  and  "V^  one  writes 

[%?*“■ r t(3V-  u)L 

it 

Ec  and  'Pi  denote  values  of  E  and  ~P  at  the  initial  time. 
Equation  (4.2  -  1)  may  be  written 

-f  2  $(i-  ^r^^V-v)+(sr-sv)j  +  Jf- 


(4.2  -  4) 


(4.2  -  5) 


By  differentiating  equation  (4.2  -  2)  with  respect  to  ^  and  adding  to 
the  resulting  equation  the  product  of  and  equation  (4.2  -  2)  one 
obtains 


(4.2  -  6) 


Substituting  equation  (4.2  -  4)  into  equation  (4.2  -  3)  and  then  differ¬ 
entiating  with  respect  to  one  finds  that 

°f'~ 

-/  V*  \f’~  1  \ 


tc6  w  ' 'Jet.  (4.2-7) 

Substitution  of  (4.2  -  6)  and  (4.2  -  7)  into  equation  (4.2  -  5)  yields 


where 


(4.2  -  8) 
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*[-  FrF^*2^- 

f.j/ffif*  (4  2-9) 

'  ?<(-'  1  ~ 

The  function  Fz  contains  temperature  and  pressure,  material  pro- 
perties,  and  the  strain  at  the  initial  time,  hence,  is  a  function  of 
time  and  ^  . 

The  first  order  partial  differential  equation  (4.2  -  8)  is  readily 


solvable  by  the  method  of  separation  of  variables. 

Employing  the  same  techniques  as  used  in  the  evaluation  of 
it  can  be  shown  that 


Thus 


"V=TT-e 


-U 


Fclir  A 


'£t. 


(4.2  -  10) 

Replacing  equation  (4.2  -  10)  for  the  value  of  "V  in  equation  (4.2  -  3), 
one  writes  that 


Ft 


where 


-fat 


(4.2  -  11) 


Ft  ,  Al—  I fl- 

Jj"f  f,  f%Jt  ,  ) 

[-{&"  1) 


(4.2  -  12) 


Using  the  method  of  separation  of  variables  on  equation  (4.2  -  11)  and 
integrating  one  obtains  v 

iyJHois  the  value  of  TT  at  time  ~t  and  at  which  is  the 

edge  of  the  s-iell,  consequently  it  is  determined  by  the  edge  conditions. 
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For  example,  if  the  edge  is  of  fixed  support,  then  = 

With  value  of  established  can  be  evaluated  from  equation 
(4.2  -  10).  The  membrane  forces  and  moments  are  subsequently  evaluated 
by  equations  (3.1.1  -  5)  through  (3.1.1  -  8). 


V.  BENDING  ANALYSIS 


5.1  General  Shell 

To  more  accurately  describe  the  effects  of  edge  conditions  the 
following  derivation  includes  the  changes  in  curvature  of  the  middle 
surface. 


5.1.1  Governing  equations 

Let  -€f-  tty  and  £©  y  =  -  ?X&  where  £y  and  €■*.  are  the 

middle  surface  strains  and  ~y.x  and  ^  are  the  changes  of  curvature. 


c  =J_fQL-uo\ 


^2  = 


-  -f-  /  ir  cot  -  u/ 1 

Jr-L  \  * 


*~k(i 


hn\ 

A 


(5. 1.1-1) 
(5. 1.1-2) 
(5. 1.1-3) 
(5. 1.1-4) 


For  the  Kelvin-Voigt  body,  the  stress-strain  relationship  is 


2 


2  c 

Proceeding  with  the  method  outlined  by  Timoshenko  where 

/A 


then 


Hj,  - rfe.$)  -  /i  (2% *x„) -  r(*ih a  fa -f) 

tip-*  A  (z^t^)  -Trfefy+Xt')  -  + TT  fa  fa) 


where 


r(-)  -  4%  e  *r'^»  +2<*0-fn.)  '(•■■) 

AN  *  He^(£)'[IP‘M¥-r"*%-]2&L*  ■frW&O"') 

7 rf  ) 

Assuming  Tm  and  Tfc  independent  of  H  then 

3/” fijr ^ ^  ^(3* (3ai^  (5. 1.1-5) 

JT  (rurfufaj  f+j  $$  *>r/]-r{3«Tj+A  (5 .  l .  1  -6) 

Eliminating  «/  from  equations  (5. 1.1-5)  and  (5. 1.1-6)  we  have 

5/-  j Jp  )-fz (2fJ6-ty)+3h 1'1_7) 

By  differentiating  equation  (5. 1.1-6)  with  respect  to  <&  and  using 
equation  (5. 1.1-7)  to  eliminate  then 

Jrfas f)~  )- 

(5. 1.1-8) 

let  -r-M«+  *%)  and  then,  from  the  equations  of  equili¬ 
brium  for  p  =  0,  ~ 

~LTcot<£ 

//. _ LJS- 

™  n  ¥ 


(5. 1.1-8) 
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With  the  above,  equation  (5. 1.1-8)  becomes  the  first  governing  equation. 


A 


HKihi-f]  S- 


(5. 1.1-9) 


To  obtain  the  second  governing  equation,  substitute  the  expressions 
for  and  into  the  equation  of  equilibrium 

■^(a  cost^ -A  ir=  ° 

then 


(5.1.1-10) 


To  eliminate  the  middle  surface  strains  from  (5.1.1-10)  consider  the 
following  equation  of  equilibrium 


I 

(5.1.1-11) 


Substituting  the  expressions  for 


and 


into  (5.1.1-11)  we  have 


(5.1.1-12) 


With  equations  (5.1.1-10)  and  (5.1.1-11)  we  now  have  the  second  govern¬ 
ing  equation 


(5.1.1-13) 
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then  the  governing  equations  become 

-*r[Lfr)+fir]+3t (5-l-l.13a) 

L(v)-y.V~=j[rfr)-  (5.1.1-9a) 


5.2  Conical  Shell 

for  the  conical  shell 


ft  =  twY 


■  -  o 


f, 

Pdd>=2t 

,,  ,  ,  vrt  2Y-)  .  rf-)  If.. A 7 

then  £.(*••/  becomes  /.(•■-)-  J t  *  J 

and  the  governing  equations  are 


-2n  [l(v)J  +  2/1  (-fr[/[i(r)}]J^V7»r 
L(v)-ir(v ) 


(5.2-1) 

(5.2-2) 


where  ~V~  is 

now 

Solving  equations  (5.2-1)  and  (5.2-2)  for  V  then 


-2JT  {u(r)  ] +2/1  {■£(/ \[u(r)flj-  ir(v) 


(5.2-3) 


Assuming  a  product  solution  Tft)  then 

yML  jriii  _ .  * 

&  -2JT(T)+2A  “  (5.2-4) 

where  the  constant  /i4  must  be  /z/t-iP)  /j*  for  the  solution  to 

conform  to  the  initial  condition.  ' 

The  solution  as  a  function  of  1  is 

M)  -  c,  [?,(?)+  f 
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where 


The  equation  resulting  from  (5.2-3)  describing  the  solution  as  a 
function  of  time  is  a  first  order  differential-integral  equation  which, 
when  differentiated,  yields 

[i+4-4]  $  +[[(3^AfAA?i(-pu)et+2AMfl!- $+ 

-> [ ('?**)¥“+  [fMA  - 

+A$  (5.2-5) 

where 


Az  -  V  (A)z  [(*£]+  *&-] 

A  is  Go^Tt 

3  ZfySutH 

a+  -  a<»  ^  -  “w  -^y 

4  -  {[ m+2]  «  Ay-  ¥■} 


By  expressing  the  coefficients  of  T  and  its  derivatives  as  Taylor  series 
then  T(t)  may  be  expressed  as  a  power  series 


77*)  -  «. 


<2,^  -f  ... 
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where 


<k -  («,&-m]+‘ilwxs*])/w  s  / 

4j  =  - 


a  = 

«n 


.«£^v*4SF.  «v,4f>,» 


-h 


■■*  X-fe-fXti]} 


*1 


(•-•)  % 


4  =  (l-A*)Ai'AAr 

R  =  \4-A/-Ai_ 

Bz  =  2^,^  (3-A)ArAtAr 

B4  =  AtAr-\Ai 

Bs'iti-AMr 


The  first  two  coefficients  («e, 4,  )  are  then  determined  from  the 
particular  boundary  conditions.  '  / 

5.2.1  Simplified  Conical  Shell 

Since  the  programmed  temperature  rates  in  the  test  schedule  pro¬ 
duced  a  temperature  difference  across  the  thickness  of  less  than  one 
degree,  the  following  analysis  neglects  the  differential  (Tt  =  0). 

The  operators  defined  in  section  5.1.1  become 

rn  -  4  } s 

A  o 

Tf(-)  =  zjie**  Zg-+4.0-fiK)(--) 
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,  then  equation  (5.2-3) 


(5. 2.1-1) 


Let  rH  =  Z\%£  6c  (-&)(•  •’ )  ]  ,  then  equation  (5.2-3) 

becomes 

Xr[aM]^VrM 

then,  with  a  product  solution,  we  have 

rM[u(r)+f,$}-*  (5. 2.1-2) 

where  is  the  same  as  that  in  section  5.2  and  T(t)  must  be  deter¬ 
mined  from  the  boundary  conditions. 

5.2.2  Boundary  Conditions 

Consider  the  edge  of  the  shell  fixed  against  rotation  and  expan¬ 
sion,  then  at  Q=f0 

Fom  all  i-  9  9_T\ 


iO  Fbtt  ALL  f 


(5. 2. 2-1) 


(5. 2. 2-2) 


Then 


V.r.-° -13.*  t 


3Sr(^2He-f/r2sr(«Tjl 


(5. 2. 2-3) 
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where 


//ru«y[jt  ffi]  r(rj- 

y,  -  2iuTr  Ta  _  ip  +  4p->A„  ff*2fr'  ffijr 'frj -  fi"»* 


Applying  equation  5. 2. 2-2  then 

r(rh  -+ 


*k. 

>2 


.2&: 
i* 


pt.r*,r-ur[%6 

£-${*rY  [  ?t«  -  f  ,»  4  ?/»*  jZ  g'»  ,  *g/«»  s 

'I*.  Fj^  ( J7fe7(  >*jrirJ] 

The  inner  surface  stresses  are  then  written  as 

7’^Z'  f1  ■  r/i  - 

-  tar  y[f-  f,  - 

r,-f^,3arr[f-  f+*f->4{  ff*  f)]r(rh 

-fj*»r-3<*r[ f*  &+A,(jLf2p.  4*)]r(r) 

(5.2. 2-4) 

Consider,  now,  thehinged  expansion  edge;  that  is,  an  edge  free  to 
rotate  completely  and  to  expand  to  some  degree. 


Then,  at  /=■  (0 


Foa.  ALL  ■& 


6Z  -  <*eT^ 


(5. 2. 2-5) 
(5. 2. 2-6) 


where  crtf7ff  denotes  the  expansion  of  the  mount.  So  with 

«,  -  2-fj]  fr(r) 

then  with  (5. 2. 2-6)  ^ 


L - 1 


i,. 

J2L. 


23f- 

jT 


K  _  3.  .  2zL 

*  x  ;; 
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Applying  equations  (5. 2. 2-3)  and  (5. 2. 2-6)  then 

rr/rn)  ^  wf-CS [%C 


r(r ')■ 


5.3  Spherical  Shell 

H 

Neglecting,  again,  the  temperature  differential,  equations 
(5. 1.1 -9 a)  and  (5. 1.1 -13a)  become,  for  the  spherical  shell  of  radius  a, 


where 


Lfp)-£iv-~-$srM 

-4r{L(v)+-h.-rj.-rr 


(5.3-1) 


(5.3-2) 


The  exact  solution  of  equations  (5.3-1)  and  (5.3-2)  is  a  form  of 
the  hypergeometric al  series  whose  convergence  is  unsatisfactory  for 
{&/£  >  to  ,  The  following  solution  will  be  an  approximation  developed 

by  Timoshenko 2 5  for  shells  whose  angle  is  not  small. 

Let  &i  =<5^ and  ~V^=~V  then  the  first  derivatives  in 

equations  (5.3-1)  and  (5.3-2)  disappear.  Neglecting  the  quantity  in 
comparison  with  its  second  derivative  we  have 


-VfM 

?fizL-JgQ, 

Txfi2  ^ 


(5.3-3) 


(5.3-4) 


the  solution  of  which  is 


Qf.  -  e~* 


(5.3-5) 


where  'f  is  the  shell  half  angle  (not  small) ,  /}  -  -Jyj  and  Y 


is  an  arbitrary  constant. 
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I 


I 

| 

Then 

=  2 ar( r-ffl -  4 

tf,  -  4 

=  $*%>£+)]  2[a*(w+r)+&(’ty+)ri]TM 
/%  -  4x^\f/tJM-f)] 2[ll+2£orfr-<fl}t* (zfrrj+swfiMjTfcl 


Employing  the  fixed  edge  conditions  on  a  hemisphere  the  inner 
surface  stresses  are  then 


<y 


-v)/*  /  1  —nr 

-Mi1 /W**  ¥)  -/2-^Vf -r)}*4t+ij\TM  - 

vp^-t) L  «  J  , 

[[,f  hr(-T-fj]»(^  ? 


where 


2//) 

Tk)  -  -  sr£»  -  <-//&*  eT'2t^hr^T-l) 


With  the  hinged  expansion  edge,  the  inner  surface  stresses  are 


$^S/n(X-jr) 

% ’vi/hZ/f-r)  l2‘x ar(f ■& + 

-  in  1  {Hy "frtifm 


where 


I  •  — 

4?/tr'£(«z-«rey  )  1 

'  r  ✓  j 
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5.4  The  Simplified  Conical  Shell  as  a  Maxwell  Body 
For  a  Maxwell  body 

-*%[&*+$- *-r~] 

Proceeding  as  in  section  5.1  and  with 


tk  =  a  S«rt  $ -If*  A‘(~  f'  f-j]  rizfj- 

d, = uary[f-  2-p-+ +4j-  f'-  ^ 

ying  the  fixed  edge  conditions,  then  r<f,(l-AT*)  , 

*  mm 

■■  ' r'"'  '^I¥¥:¥-4¥¥1W1 


o't 

cJt 


For  a  constant  temperature  rate  T =  Rt.  The  integral  then  becomes 


e**'  Jk 

T> -  At 


For  <0.  /-Jr~ 


,  ^  Let  O.i 

After  t  ?  <r 
compared  to  the  exponential  so 


the  exponential  is  essentially  unity  so 

[A‘ e  <*{'-&)&  ■,  (*< 


Let  O.I 

the  quantity  (1-  6  Rt)  can  be  considered  constant 


)UM*e ^(,-HR  fc (-££&&)-*('■')} 

*  .  ....  9m  - —  /  \ 


where 


i'- 


and  Q  (•••)  is  defined  as 
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f*K  i>i' 


The  time  function  is  then 


**  *<*' 


VI.  EXPERIMENTAL  RESULTS  AND  VERIFICATION 
OF  ANALYTICAL  METHODS 


6.1  Basic  Approach  .to  the  Test  Program 

In  the  early  stages  of  the  current  investigation,  it  was  decided 
that,  for  experimental  verification  of  theoretical  methods  of  predicting 
the  viscoelastic  response  of  shells  of  revolution  subjected  to  normal 
pressure  and  elevated  temperature  and  thermal  gradients,  records  of  the 
strain  as  a  function  of  time  would  be  necessary.  Accordingly,  steps 
were  taken  to  find  appropriate  means  of  making  such  elevated  temperature 
measurements. 

The  objectives  of  the  experimental  investigation  were  as  follows: 

(1)  To  obtain  strain  versus  time  curves  for  normal  pressures 
and  applied  temperatures  and  thermal  rates. 

(2)  To  determine  the  mode  and  progression  of  failure,  and  the 
load  condition  necessary  to  cause  failure  under  load  and 
transient  temperatures. 

(3)  To  verify  or  disprove  the  method  of  analysis. 

6.1.1  Test  conditions 

The  experimental  program  for  verification  of  a  method  of  analysis 
does  not  necessarily  duplicate  any  specific  design  conditions.  It  is 
required  that  the  tests  cover  the  parameters  included  in  the  method  of 
analysis.  These  parameters  should  extend  over  a  wide  range  of  the 
values . 

For  convenience,  a  uniform  normal  pressure  of  12  psi  is  chosen  to 
be  applied  simultaneously  with  various  combinations  of  temperature.  A 
maximum  test  temperature  of  500°F  is  chosen,  based  on  a  level  of  tem¬ 
perature  sufficiently  high  to  introduce  material  property  variation 
and  buckling. 
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The  tests  are  designed  to  eliminate  as  many  extraneous  variables  as 
possible.  However,  in  the  present  study  it  is  impossible  to  eliminate 
parameters  such  as  edge  constraints,  and  discontinuities  of  temperature 
due  to  mass.  These  variables  are  difficult  to  evaluate.  Xn  the  corre¬ 
lation,  it  is  not  clear  whether  the  method  of  analysis  is  at  fault  or 
the  accuracy  of  the  experimental  data  is  questionable.  The  interrelation¬ 
ship  of  the  above  mentioned  inseparable  parameters  and  their  influence 
on  the  overall  stress  distributions  are  not  satisfactorily  established. 

6.1.2  Number  of  Tests 


The  number  of  tests  is  an  important  factor  in  the  design  and  cost 
of  an  experimental  investigation.  Usually  repetition  of  test  conditions 
with  several  specimens  is  made  to  include  the  effect  of  test  scatter. 

A  more  economical  approach  is  to  test  each  specimen  under  different 
conditions  and  to  compare  the  test  results  with  theoretical  results  over 
the  range  of  the  variables.  This  approach  permits  checking  over  a  large 
range  of  values  with  a  minimum  number  of  tests  and  test  specimens. 

The  approach  employed  in  this  investigation  is  designed  after  the 
latter  method;  however,  where  certain  doubts  existed,  the  test  conditions 
are  repeated  with  another  model. 

6.1.3  Test  Specimen 

The  size  of  the  test  specimen  is  another  important  consideration 
particularly  in  elevated  temperature  structural  research.  In  such 
investigations,  large  specimens  are  required  to  eliminate  certain  test¬ 
ing  difficulties;  for  example,  temperature  gradients  are  established  by 
temperature  difference  along  the  length  of  a  specimen.  If  the  specimen 
were  too  small,  temperature-sensing  elements  would  be  too  closely  spaced 
along  the  specimen  with  probable  disturbance  of  the  required  gradient. 

The  design  of  the  specimen  is  based  on  the  following  specifications: 

(1)  Geometry  -  conical  and  hemispherical  shells  are  used  analyti¬ 

cal  and  testing  simplicity. 

(2)  Radius  -  a  base  radius  of  12  inches  is  chosen  as  a  con¬ 

venient  specimen  size. 

(3)  Thickness-  the  thickness  is  based  on  a  radius-to-thickness 

ratio  corresponding  to  a  convenient  anticipated 
buckling  stress. 

(4)  Material  -  Aluminum  is  chosen  so  that  material  property  varia¬ 

tion  with  temperature,  creep,  and  buckling  occur 
within  a  convenient  temperature  and  load  range. 

This  is  of  primary  importance  in  the  design  of  the 
test  facility  and  the  instrumentation. 
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6.2  Test  Procedure 


Prior  to  an  actual  test  run,  an  unrestrained  heating  of  the  model 
provided  data  necessary  for  temperature  compensation  of  the  strain 
gages.  The  difference  between  these  data  and  the  actual  test  data  com¬ 
pensated  the  effects  of  a  varying  coefficient  of  linear  thermal  expansion, 
coefficient  of  resistivity  of  the  gages,  and  coefficients  of  linear 
thermal  expansion  of  the  gage  cement  and  gage  material.  The  compen¬ 
sated  strain  data  is  then,  in  effect,  the  total  strain  minus  the  linear 
thermal  expansion. 

With  the  strain-temperature  data  recorded  and  the  model  cooled, 
the  shell  is  bolted  to  the  steel  mounting  plate.  The  selected  time- 
temperature  function  is  set  on  the  temperature  programmer. 

Instead  of  using  an  ice-water  bath  for  the  cold  junction  of  each 
thermocouple,  a  separate  thermocouple  employing  an  ice-water  bath 
provides  the  temperature  of  the  recording  thermocouple  input  to  the 
recorder.  This  temperature,  which  is  monitored  on  a  precision  millivolt 
potentiometer,  is  the  effective  cold  junction  temperature  of  the 
recording  thermocouples.  The  strain  gage  bridges  are  balanced  and  zero 
points  of  all  instrumentation  and  reference  tanperatures  of  thermo¬ 
couples  are  noted.  The  oven  is  placed  over  the  model  and  the  power  con¬ 
troller  is  connected  to  the  lamps.  The  model  is  evacuated  to  conform 
to  the  prescribed  external  load,  after  which  the  recording  instruments 
and  the  temperature  programmer  are  started. 

After  completion  of  the  transient  portion  of  the  test  program, 
the  temperature  programmer  is  placed  in  a  "hold"  condition  to  maintain 
the  steady-state  temperature  for  the  desired  period  of  time.  Usually 
the  steady-state  is  maintained  for  five  minutes  except  for  collapse 
studies  in  which  case  thirty  minutes  are  provided. 

When  the  test  is  completed,  the  power  controller  is  shut  off,  the 
model  is  allowed  to  cool  and  then  unbolted  for  the  next  unrestrained 
run. 


6.3  Analysis  of  Test  Results 

The  experimental  results  are  presented  in  Fig.  6.5-1  through  6.5-18 
and  included  in  Appendix  C. 

The  low  thermal  resistance  of  models  requires  a  low  gain  condition 
on  the  power  controller  to  prevent  a  pulsing  heat  input.  As  a  result 
of  the  low  gain,  the  steady-state  temperatures  are  40°F  to  60°F  lower 
than  that  programmed  in  the  initial  tests. 

The  base  temperatures  are  considerably  lower  than  that  of  the  rest 
of  the  shell  due  to  heat  conduction  to  the  steel  mounting  plate.  In  an 
effort  to  remedy  the  low  temperature  at  the  base,  the  feedback  thermo¬ 
couple  for  Zone  3  was  moved  closer  to  the  base.  This  produces  a  shell 
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temperature  40°F  to  80°F  higher  than  the  program  temperature,  but  the 
temperatures  near  the  base  are  from  100°F  to  150°F  lower  than  the 
steady-state  temperature.  However,  the  resulting  temperature  gradient 
occupies  only  the  lower  three  inches  of  the  shell. 

The  data  from  the  differential  thermocouples  are  merely  qualita¬ 
tive.  Calculations  show  a  thermal  gradient  across  the  shell  thickness 
of  1°F  for  a  programmed  temperature  rate  of  20°F/sec  whereas  the  data 
indicates  30°F.  The  large  experimental  values  are  due  to  a  poor  thermal 
bond  between  the  insultated  junction  and  the  surface  of  the  shell, 

A  characteristic  feature  of  all  tests  is  the  relaxation  of  stress 
as  steady-state  temperature  conditions  are  maintained.  The  reduction  of 
compressive  stress  is  due  to  the  expansion  of  the  steel  mounting  plate 
and  a  relaxation  of  the  material .  The  relative  contribution  of  each 
condition  depends  upon  the  temperature  of  the  base.  To  separate  the  two 
effects  requires  insulation  of  the  base.  However,  it  is  felt  that  the 
major  factor  in  the  reduction  of  compressive  stress  is  that  of  material 
relaxation.  This  conclusion  is  supported  by  some  of  the  zonal  tests. 

In  those  tests  where  the  base  temperature  was  relatively  low,  say  about 
120°F,  in  which  case  yielding  of  the  base  should  not  be  a  problem,  the 
compressive  stresses  along  the  model  slant  height  decayed  although  the 
compressive  stresses  near  the  base  increased  with  an  increasing  steady- 
state  temperature. 

As  is  expected  the  stresses  are  compressive  throughout  the  shell 
with  those  in  the  circumferential  direction  being  of  the  greater  mag¬ 
nitude.  Maximum  values  are  observed  at  the  base,  and  the  stress  decreases 
rapidly  along  the  meridian  fbr  the  lower  section.  The  upper  regions  of 
the  shells  experience  approximately  the  same  stress  magnitude.  The 
meridional  stresses  are  similar  to  the  circumferential  stresses  but  with 
less  pronounced  maxima  and  minima. 

Collapse  was  encountered  in  a  1/16"  thick  conical  shell  during  an 
applied  heating  rate  of  20°F/sec  (to  500°f)  with  a  surface  pressure 
equivalent  to  22"Hg  differential  pressure.  The  same  model  had  experienced 
the  same  program  previously.  The  data  are  presented  in  Table  C-ll. 
Approximately  at  80%  of  the  transient  period  gage  4M  showed  an  immediate 
change  in  strain.  The  change  is  believed  to  be  caused  by  a  slight 
buckle  at  the  support.  No  further  buckling  occurred  until  the  transient 
condition  had  been  sustained  for  five  minutes  when  all  gages  registered 
an  immeuiate  jump.  Thereafter  gages  3M  and  3P  showed  increasing  com¬ 
pression  to  five  and  one-half  minutes  when  complete  collapse  occurred. 

The  collapsed  cone  has  five  non  axisymmetric  nodes  with  the  wall 
in  the  region  of  the  base  almost  flat  against  the  support.  After  col¬ 
lapse  the  support  for  instrumentation  wiring  interf erred  with  progression 
causing  a  tilt  of  the  apex  toward  the  more  seriously  buckled  area. 

See  Fig.  6.3-1. 
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6.4  Simplified  Strain-Displacement  Method 

Compared  in  Fig.  7.1-4  and  7.1-5  are  the  experimental  data  for  a 
1/16"  thick  cone  for  20°F/sec  to  500°F  with  11.5  psi  pressure  and  the 
theoretical  stress  calculated  from  the  equations  presented  in  section 
3.1  for  a  similar  program.  As  seen  from  the  figures  the  theoretical  stress 
behavior  is  opposite  to  that  of  the  experimental.  This  behavior  is 
believed  due  to  neglecting  the  change  in  curvature  in  the  strain  ex¬ 
pressions  thereby  relegating  the  analysis  to  a  non-bending  case  (mem¬ 
brane  theory) . 

6.5  Bending  Analysis  Neglecting  the  Temperature  Differential. 

Compared  in  Fig.  6.5-1  through  Fig.  6.5-18  are  the  theoretical  and 
experimental  data  for  one  test  on  a  1/32"  thick  cone,  five  test  on  a 
1/16"  cone,  one  test  on  a  1/8"  cone:,  and  two  tests  on  a  1/16"  thick 
hemisphere.  All  theoretical  data  are  calculated  for  the  hinged  expan¬ 
sion  edge  boundary  condition,  i.e.,  the  mounting  plate  is  considered  to 
assume  the  temperature  of  the  edge  of  the  model  and  the  edge  offers  no 
resistance  to  rotation  in  the  merional  direction.  A  characteristic  of 
all  theoretical  data  is  a  large  stress  value  during  the  early  time 
periods  when  the  response  should  be  nearly  elastic.  The  time  periods 
are  1/3,  1/2,  and  end  of  transient. 

The  Kelvin-Voigt  body  may  be  pictured  as  a  spring  and  a  dashpot 
in  parallel  with  each  element  having  the  same  displacement  and  the  re¬ 
sulting  force  being  the  sum  of  that  in  each.  Since  the  applied  boundary 
condition  of  a  restrained  edge  is  essentially  an  applied  strain  rate, 
the  large  stresses  are  developed  in  the  viscous  member  conforming  to 
the  boundary  condition.  At  such  time  when  the  temperature  is  increased 
to  cause  an  appreciable  effect  on  the  coefficient  of  viscosity  (the 
coefficient  varies  inversely  as  the  temperature)  the  viscous  resistance 
becomes  negligible  and  the  response  is  one  which  is  mainly  elastic.  ,  % 

From  section  5.2  the  thermal  stress  is  governed  by  the  term  % tKT~ 
so  with7„'v/0tt  and^j/v/a*  the  viscous  stress  is  necessarily  large  1?‘ 

for  low  temperatures. 

At  higher  temperatures  the  behavior  of  the  theoretical  curves  are 
similar  to  that  of  the  experimental,  that  is,  an  increasing  compressive 
stress  toward  the  edte  increasing  iso  with  time. 

Since  the  stresses  are  calculated  at  the  inner  surface  they  show 
the  effect  of  the  bending  moment  close  to  the  edge.  The  theoretical 
bending  moment  accounts  for  the  majority  of  the  stress  in  the  meridion¬ 
al  direction.  The  moment  changes  rapidly  with  distance  from  the  edge 
and  depends  considerably  upon  the  boundary  conditions. 

The  experimental  boundary  conditions  cannot  be  assumed  ideally 
hinged  or  fixed.  Due  to  heat  conducted  from  the  edge  to  the  mounting 
plate  the  edge  will  expand  but  not  to  the  extent  assumed  in  the 
expansion  boundary.  Also  the  edge  may  rotate  but  is  not  ideally  hinged. 
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Actual  edge  conditions  would  require  an  analysis  investigating  a 
flanged  edge. 

The  comparison  in  Fig.  7.1-4  and  7.1-5  shows  the  effect  of  boundary 
conditions  on  the  inner  surface  stress.  The  two  different  sets  of 
boundary  conditions  should  bracket  the  actual  situation.  In  the  fixed 
edge,  the  meridian  may  not  rotate  and  the  edge  may  not  expand.  For  the 
hinged  expansion  edge,  the  meridional  moment  at  the  base  vanishes  and 
the  ed^e  may  expand  with  the  steel  plate  (  ^aluminum  ,  ^steel 

rj  ~)uo  ).  The  fixed  case  produces  tensile  stress  and  the  hinged  case 
compressive  stress  near  the  edge  on  the  inner  surface  of  the  shell  in 
the  meridional  direction.  The  effect  on  the  middle  surface  is  shown  in 
the  following  figure. 


The  circumferential  stress  depends  mainly  upon  the  degree  of 
restraint  at  the  edge  and  is  proportional  to  the  difference  between 
the  expansion  of  the  shell  and  the  expansion  of  the  mount. 

The  experimental  values  compare  more  favorably  with  the  theoretical 
stress  for  the  hinged  expansion  edge.  The  maximum  experimental  circum¬ 
ferential  stress  is  about  607»  of  that  for  the  hinged  edge  while  40%  of 
that  for  the  fixed  edge.  The  maximum  experimental  meridional  stress  is 
about  75%  of  that  for  the  hinged  edge  while  the  fixed  edge  stress  is  of 
the  same  order  of  magnitude  but  of  opposite  sign.  While  the  experimental 
stresses  should  fall  between  those  of  the  two  boundary  conditions  it  is 
seen  that  the  measured  stress  is  still  lower  than  either.  The  Kelvin- 
Voigt  body  does  not  consider  the  relaxation  of  stress  at  constant  strain 
at  high  temperatures  or  the  yielding  of  the  material  for  high  stress. 
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6.6  Comparison  of  the  Maxwell  and  Kelvin -Voigt  Bodies 

In  Fig.  6.6-1  and  6.6-2  the  viscoelastic  response  of  a  1/16"  thick 
cone  subjected  to  20°F/sec  to  500°F  is  compared  for  a  Maxwell  body  and 
a  Kelvin-Voigt  body.  The  early  response  of  the  Kelvin-Voigt  body  is, 
again,  of  a  highly  viscous  nature  while  the  Maxwell  body  response  is 
essentially  elastic.  At  higher  temperatures  when  the  Kelvin-Voigt  body 
becomes  elastic  the  Maxwell  response  is  fluidic  due  to  the  applied 
strain  rate  being  absorbed  in  the  viscous  component.  The  theoretical 
stress  is  essentially  the  product  of  a  function  of  slant  length  and  a 
function  of  time.  With  the  function  of  slant  length  being  the  same 
for  both  bodies  the  time  function  is  presented  in  the  following  figure: 


T(t) 

Kelvin 


Moreover,  at  the  steady  state  condition  the  Kelvin-Voigt  body  maintains 
a  constant  stress  while  the  Maxwell  body  will  decay  exponentially. 

Compared  in  Fig.  6.6-3  and  6.6-4  are  the  results  of  the  temperature 
dependent  properties  analysis,  the  temperature  independent  properties 
analysis,  and  the  temperature  dependent  properties  elastic  analysis. 

The  property  values  for  the  temperature  independent  analysis  were  chosen 
at  a  temperature  corresponding  to  10  seconds.  The  resj>onse  of  the 
temperature  independent  analysis  depends  upon  the  choice  of  the  coeffi¬ 
cient  of  viscosity;  if  the  temperature  is  low  the  response  will  be 
viscous  and  if  the  temperature  is  sufficiently  high  the  response  will  be 
elastic. 

The  variation  of  the  elastic  solution  from  that  of  the  viscoelastic 
solution  at  high  temperatures  is  due  to  the  choice  of  Poisson's  ratio 
being  1/3  for  the  elastic  analysis  and  \  for  the  viscoelastic.  Again 
the  experimental  stress  is  lower  than  any  theoretical  prediction  which 
is  believed  due  to  the  obscure  boundary  conditions  and  a  yielding  at 
the  edge. 
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VII. 


SUMMARY  AND  CONCLUSIONS 


The  experimental  results  in  some  instances  are  in  disagreement 
with  the  temperature  dependent  viscoelastic  analysis.  This  is  gen¬ 
erally  true  for  the  early  transient  stage  and  the  steady-state  stress 
conditions.  In  an  effort  to  explain  these  differences,  a  number  of 
questionable  elements  of  both  the  theoretical  analysis  and  the  exper¬ 
imental  investigation  were  studied.  These  points  are  discussed  in 
the  following  paragraphs,  and  when  possible,  recommendations  to  eli¬ 
minate  the  individual  effects  are  given. 

7.1  Theoretical  Analysis 

7.1.1  Viscoelastic  models 

Employed  in  the  analysis  is  the  assumption  of  the  material's 
response  as  that  of  a  simple  Kelvin -Voigt  body.  A  much  better 
approximation  would  be  the  development  of  a  model  consisting  of  a  number 
of  Kelvin-Voigt  elements  in  series.  The  difficulties  in  the  analysis 
during  the  early  transient  and  the  steady- state  condition  would  benefit 
from  a  more  descriptive  model. 

From  the  comparisons  in  Fig.  6.6-3  and  6.6-4  it  appears  that  the 
elastic  analysis  would  sufficiently  describe  the  stress  response  in 
the  models  tested  as  compared  to  the  viscoelastic  analysis  presented. 
However,  neither  analysis  predicts  the  steady- state  response. 

7.1.2  Effect  of  material  properties 

In  the  present  investigation,  an  empirical  expression  for  was 
established  by  employing  some  limited  data  established  by  ALCOA  at 
higher  temperature  ranges  for  aluminum  alloys.  Because  the  viscosity/ 
temperature  relationship  is  based  on  such  limited  data,  three  values 
were  used  in  the  calculations.  The  effect  of  varying  ^  is  demonstrated 
in  Fig.  7.1-1  and  7.1-2.  Changes  of  an  order  of  magnitude  are  employed. 
As  seen  from  the  figures,  the  only  effect  is  to  advance  or  retard  the 
transition  from  a  viscous  response  to  an  elastic  response. 

Although  some  error  may  exist  in  the  choice  of  G  and  ^  ,  it  is 
felt  that  both  are  minor  and  have  little  effect  upon  Che  results. 

RECOMMENDATION:  More  reliable  data  on  the  dependance  of  thermal 

properties  of  alloys  with  temperature  are  required  for  accurate 

analyses.  In  particular,  viscosity/temperature  data  on  alloys 

below  the  melting  point  should  be  established. 
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7.1»3  Inaccuracies  in  the  applied  boundary  conditions 

Agreement  of  the  boundary  conditions  employed  in  the  analysis  with 
those  acutally  produced  in  a  test  program  is  difficult  to  achieve. 
Although  such  disagreement  is  present,  it  is  felt  that  the  effect  is 
localized  to  the  edge  and  apex  (in  the  care  of  cones)  regions  as  far 
as  the  stresses  in  the  shell  are  concerned. 

7.2  Experimental  Study 

Aside  from  the  usual  sources  of  errors  encountered  in  experimental 
stress  investigations  and  data  reduction,  the  investigators  feel  that 
the  following  items  are  major  possible  contributors  to  the  existing 
disagreement  between  theoretical  and  experimental  results. 

7.2.1  Inaccuracy  in  correcting  strain  data  to  eliminate  temperature 
effects  in  strain  gages 

There  is  much  to  be  done  to  perfect  high  temperature  strain  mea¬ 
suring  techniques.  The  range  of  temperatures  of  interest  is  subject 
to  a  number  of  technical  problems  both  in  terms  of  model  preparation 
and  testing  procedures  for  T>400°F.  In  fact,  the  experience  gained  in 
this  investigation  indicates  that  it  requires  six  to  eight  months  to 
develop  techniques  for  such  an  investigation. 

RECOMMENDATION:  High  temperature  strain  gages  which  require 
less  complicated  application  procedures  are  needed.  These  strain 
gages  should  be  temperature  compensating.  Since  the  completion 
of  the  present  investigation,  other  high temperature  strain  gages 
have  been  brought  to  the  attention  of  the  investigators.  It  would 
be  of  interest  to  note  whether  these  newer  developments  eliminate 
the  difficulties  of  application,  curing,  and  measuring  as  ex¬ 
perienced  in  this  investigation. 

7.2.2  Effect  of  model  support  on  temperature  distribution 

The  steel  plate  which  supported  the  oven  and  the  model  had  consider¬ 
able  effect  on  the  model  temperature  near  the  support.  Although  steps 
were  taken  to  elirainate  the  plate  as  a  heat  sink,  its  effects  are  evident 
in  the  experimental  results.  The  existence  of  such  a  heat  sink  is  per¬ 
haps  the  major  difficulty  of  the  equipment.  It  created  an  undesirable 
thermal  gradient  along  the  shell  length.  Otherwise  the  equipment  per¬ 
formed  well.  The  equipment  was  not  employed  to  the  maximum  capacity. 

RECOMMENDATION:  For  more  effective  use  of  the  present  equipment 
design,  a  ceramic  coating  should  be  used  at  the  shell-support 
connection. 
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7.2.3  Imperfection  of  the  Model 


In  Appendix  A. 3,  the  effect  of  geometrical  irregularies  of  the 
initial  shape  are  discussed. 

The  1/32"  thick  conical  shells  were  more  likely  to  contain  initial 
imperfections.  This  probably  accounts  for  the  poorer  correlation 
between  analytical  and  experimental  results.  However  for  the  thicker 
models,  not  only  was  better  correlation  achieved  but  in  several  cases 
there  was  good  agreement  of  analytical  and  experimental  results. 

On  the  basis~of  these  qualitative  results,  the  following  recom¬ 
mendation  :ls  made: 

RECOMMENDATION:  Study  should  be  made  of  surfaces  of  revolution 

which  have  unsymmetrical  initial  imperfections. 

7.3  General  Comment 

In  the  foregoing  discussion,  several  important  aspects  of  the 
investigation  are  outlined.  In  connection  therewith,  the  investiga¬ 
tors  feel  that  the  results  are  a  step  in  the  solution  of  the  complicated 
viscoelastic  shell  problem.  Nevertheless,  in  addition,  several  other 
important  aspects  of  the  analysis  bear  further  investigation;  these 
include : 


1.  Extention  of  the  theory  to  include  the  effect  of  internal 
ring-  and  meridional -frames,  and  other  surface  constraints. 

2.  Extention  of  the  method  to  include  viscoelastic  and/or 
elastic  layered  systems  of  two  and  three  layers.  The  use 
of  outer  liners  should  be  investigated. 

3.  Future  investigations  might  also  include  both  unsymmetrical 
body  forces  and  temperature  distributions. 

It  is  needless  to  say  that  each  of  these  endeavors  are  entire 
subjects  within  themselves. 

In  connection  with  Item  2,  it  is  felt  that  the  present  method  is 
too  complicated  to  be  applied  to  layered  systems.  For  this  reason, 
the  investigators  studied  the  application  of  energy  methods  to  the 
solution  of  viscoelastic  shell  problems. 
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Fig.  2.1.3- 1 

PROPERTY  VARIATION  WITH  TEMPERATURE 
FOR  PURE  ALUMINUM  AND  6061  ALLOY 


PROPERTY  VARIATION  WITH  TEMPERATURE  FOR 
PURE  ALUMINUM  AND  6061  ALLOY 

Table  2. 1.3-1 


T  (°C) 

G  x  10  ^  (psi) 

pure  aluminum*  6061  alloy** 

(poise) 

pure  aluminum*  6061  alloy 

25 

3.48 

3.80 

2. 2xl016 

8.15xl017 

100 

3.34 

3.78 

2.1X1011 

4.22xl012 

200 

3.04 

3.44 

1.3xl07 

1 . 45xl08 

285 

2.90 

2.85 

4.8xl04 

3.66xl05 

350 

2.76 

2.10 

2. 0x10  3 

9. 30x10 3 

450 

2.61 

0.65# 

43.0 

52.9 

550 

2.47 

0.35# 

2.2 

1.38 

660 

2.18 

0.20# 

0.18 

0.06** 

670 

2.18 

— 

0.14 

— 

700 

-- 

-- 

-- 

0.055** 

* 

Estimated  by  Ke 
of  the  Viscous 

.  Refer  to  Ke,  Ting-sui:  Exp 
Behavior  of  Grain  Boundaries  i 

erimental  Evidence 
n  Metal,  Physical 

Review,  Vol.  71,  pp.  533-546,  1947. 
**  ALCOA  values. 

#  Estimated  from  G-T  curve. 
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Fig.  2. 1.3-2 

MATERIAL  PROPERTY  VARIATION  WITH  TEMPERATURE  FOR  6061  ALLOY 
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FIG.  2.6-1 


MODULUS  OF  ELASTICITY  (p.s.i.  x  10  ) ,  REF.  4. 


TEMPERATURE  VARIATION  OF  MODULUS  OF  ELASTICITY 
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Fig.  6.6-3 
COMPARISON  OF  STRESSES  BASED  ON 
TEMPERATURE  DEPENDENTAND  INDE- 
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TEMPERATURE  DEPENDENT  ELASTIC 
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PBWWRII 


The  strain-displacement  relationship  for  the  general  case  are 


As  applied  to  the  conical  shell 

A-  «  f 
fz  =  V  ta*Y 


The  strain-displacement  relation  becomes 

^r’f^'wrHr 


Introduction  of  the  dimensionless  variables 

/  ^ 


gives  for  the  strains 


Since  for  a  thin  shell 


and 


it  is  obvious  that  the  terms  having  the  order  of  magnitude 
be  neglected.  Therefore, 


may 


or 


iT-  uj  (er  V 

T~ 
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Appendix  A. 2  Elastic  Analysis  of  a  Thin  Conical  Shell  of  Constant 
Thickness  under  Constant  Pressure  and  Thermal  Loads"! 


Equations  of  Equilibrium 

iu(tyt) 

( fy + &}  ta*  ir=  ^ 


Temperature  Profile  Across  the  Shell  Thickness 


Stress-Strain  Relation 


ft  *  [■*=?*+  *  %  - 


Working  Variables 


Xc/lL) 

m  -JT 
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Governing  Equations 


where  2^  =  2-,  {,  f)  and  S9  =  *9(  f  )  are  ber  and  -bei  functions  of  y  respec¬ 
tively,  and  yz  *=  -JiifX-p/fi-yi)  2  .  The  values  of  and 

A2  are  to  be  determined  by  the  edge  conditions  of  the  shell. 


Appendix  A. 3  Geometrical  Imperfections 


In  recent  years,  considerable  discussions  have  been  made  to  ex¬ 
plain  discrepancies  between  shell  theories  and  experimental  results. 

Many  investigators  feel  that  these  discrepancies  are  due  to  the  geo¬ 
metrical  imperfections  that  exist  in  the  models. 

The  present  investigators  made  a  preliminary  study  of  the  effect 
of  geometrical  imperfections  on  the  behavior  of  spherical  shells  under 
normal  pressure  to  obtain  some  understanding  of  their  effects.  The 
analysis  is  based  on  the  axi symmetrical  type  imperfections. 

As  in  the  case  of  other  investigations  by  Chen  (6),  Kaplan  and 
Fung  (15),  Budiansky  (5),  it  is  found  that  axi symmetrical  analysis 
could  not  account  for  the  total  difference  between  theory  and  exper¬ 
iment.  The.  mode  of  collapse  of  the  models  in  the  present  investigation 
gives  evidence  that  the  initial  imperfections  of  the  models  were 
probably  unsymmetrical.  Experimental  evidence  obtained  by  Kaplan  and 
Fung  leads  on  to  the  conclusion  that  an  investigation  of  the  effect  of 
geometrical  imperfections  should  be  based  on  unsymmetrical  rather  than 
axi symmetrical  analysis. 

However,  for  such  an  analysis  to  be  meaningful,  it  must  be  an 
exhaustive  study  supported  by  a  number  of  experimental  tests.  The 
study  need  not  involve  temperature.  Such  an  investigation  was  not 
attempted  as  part  of  the  present  program. 

Previous  Investigations: 

To  the  knowledge  of  the  writers,  only  five  studies  have  been  made 
about  the  effect  of  initial  imperfection  on  the  behavior  of  thin  shells. 
Four  of  the  studies  were  concerned  with  spherical  shells  and  one  dealt 
with  general  thin  shells.  In  chronological  order  the  studies  and  their 
general  results  are  as  follows: 

Mushtari  presented  the  equilibrium  equations 
of  the  middle  surface  of  a  thin  shell  which  have 
irregularities  of  the  order  of  the  shell  thick¬ 
ness.  The  relations  were  applied  to  thin  shells 
subjected  to  normal  pressure. 

Klein  used  two  parameters  in  connection  with 
experimental  data  to  provide  a  better  correlation 
of  analytical  and  experimental  results.  He 
attempted  to  show  that  by  proper  choice  of  the 
allowable  collapse  pressure,  the  effect  of  initial 
imperfections  of  the  shell  could  be  taken  into 
account. 

Gerard  and  Becker  using  an  empirical  approach 
based  on  an  "unevenness  factor"  proposed  that  the 
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geometrical  portion  of  the  initial  imperfections 
governs  the  average  behavior  of  the  spherical 
plates  which  the  residual  stress  and  6ther  fac¬ 
tors  may  contribute  to  the  scatter  in  the  experi¬ 
mental  data. 

Chen  considering  the  ax i symmetric  buckling  under 
uniform  pressure  of  a  shallow  portion  of  a  spher¬ 
ical  shell,  carried  out  a  Rayleigh-R.it z  solution 
of  the  variational  equation  equivalent  to  the 
governing  equations.  He  concluded  that  the  buck¬ 
ling  pressures  may  be  appreciably  affected  by 
the  presence  of  imperfection  of  the  middle  sur¬ 
face  and  that  its  effect  depend  not  only  on  the 
magnitude  but  also  on  the  location  of  the  maximum 
imperfection  and  the  mode  of  the  imperfection. 

Budiansky ,  treating  the  same  problem  as  Chen, 
based  his  analysis  on  the  integral -equation  for¬ 
mulation.  Budiansky  found  that  P appears  to 
approach  1  in  an  oscillatory  fashion  as  A 
(geometrical  parameter  )  increases.  Comparing  the 
variation  buckling  pressures  of  initially 
perfect  and  imperfect  clamped  shallow  spherical 
shells,  he  found  that  both  followed  the  same 
trend. 

Comparison  of  these  investigations  show  very  little  agreement  of 
the  results.  Within  a  certain  range  of  A  ,  say  =  4  to  A  =  5, 
fairly  good  agreement  has  been  observed.  Beyond  this  narrow  range, 
the  correct  theoretical  variation  of  Pcr  with  A  is  questionable. 

Method  of  Analysis 


To  study  the  possible  causes  for  discrepancies  which  exist  between 
the  present  analytical  and  experimental  results,  a  simplified  analysis 
was  made  of  thin  spherical  shells  with  axisymmetrical  imperfections. 

For  simplification,  the  study  is  restricted  to  the  hemispherical 
shells  subjected  to  pressure  only.  Following  a  method  proposed  by 
Budiansky,  we  consider  the  buckling  of  such  shells  whose  deviated  surface 
may  be  represented  by 

^  «  H[*-  (if- 6  £(*-)] 

where  e(r)  is  the  shape  of  the  imperfection,  G  is  the  ratio  of  the 
downward  initial  displacement  at  the  center  of  the  shell  to  the  rise 
H  of  the  perfect  shell. 

For  simplicity  the  shape  of  imperfection  is  taken  as 

cN  -  ['  -  tiT] ' 
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By  the  aid  of  the  geometrical  parameter 


A 


Hi-*-)  ay 


defined  by 


and  the  non-dimensional  variables 


X- 


M 

a. 


8-  (&)fi 

X  _  f  uA-i>x)a-l  -./> 

L  xzt'  J  * 

where  l/'"  is  the  stress  function,  is  the  rotation  of  an 

element  of  the  shell,  is  the  classical  buckling  pressure  of  ^a 

complete  spherical  shell  having  the  same  radius  of  curvature  R^^/zh 
as  the  given  shallow  spherical  shell. 


The  non-dimensional  equilibrium  and  compatibility  equations  for 
an  initially  perfect  spherical  shell  are 


(x£)')  -  -r  X  $  = 

(*£')- JL  -x£>  -  --le* 


The  prime  denotes  differentiation  with  respect  to  x  and 
dimensional  pressure  parameter  p  is  defined  as 

. 

where 


P 


A  non- 


The  non-dimensional  boundary  conditions  appropriate  for  a  clamped 

shell 


e(\)  =  © 


the  first  condition  states  that  the  rotation  at  r  =  a  is  zero.  The  second 
boundary  condition  fulfills  the  requirement  that  the  horizontal  dis¬ 
placement  at  the  shell  edge  must  vanish. 

For  the  initially  imperfect  shell,  the  basic  equations  take  the 

form 
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(*$')  -  ^  +  K  $~-2j>x  t  &<£+ 

(x$0  4-e*-  ^ 


2  /' 


tvMMAer  A=2l£ 


The  boundary  conditions,  of  course,  are  still  valid. 

The  solution  of  these  equations  and  the  boundary  conditions  can 
be  represented  by  two  integral  equations 

M+  &  Alfa) 

£(*.)  -  zfaj+c/fy  -2-px 

where 

y/xU  Jpfa)*{r)</r  *  £h{k;)s(t)</! 

*(*■)  -  -l#fa)*(r)Jc+j(p6,t)s/f)Sr 

/L  ~  -h£=/t 

2  -  g{,6 

2fX  (<-*)- 0-y)  ifahifal  s^* 


6(fj)  *  ke>  *  **<t'z+Xe'x tto'f  faxif) 

=  Be  I  'f  *£>i'x  -  /Sell  -fffk/'x 

H(Xt  %)  -  Btdt BE/ 'xXE/'f  {*<X<  f) 

-  Bez'%  Ket  'x  -  Bet  'fket'x  (o^f^x) 


AtM~-  £(b£s!'ABe/  X-  Set  'a Beh'x) 

=  -£( Ha! \  Bb£x + BE)  *. A  B£i'x ) 

A^Cx)  =  -(ti-y) B,U) -■^(bbh.) Be/ x- Bed  bbz'x) 

%_U)  =  (i+vjAt6)+-&(**A8a,A+*aA**'>) 

D-  ( i+v)JJbez.  ’A)  V  ( set  U)2J  -  J  (see.  A  be/  'A  -  be/  A  be*.' A  ) 
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The  solution  of  the  basic  equations  is  obtained  by  iterative  pro¬ 
cedures  in  terms  of  matrix  operations  with  the  integrations  being  per¬ 
formed  numerically  on  the  basis  of  Simpson's  rule.  For  the  particular 
application  under  discussion,  the  following  data  are  selected  for  com¬ 
parison  of  the  buckling  stress  of  a  perfect  spherical  shell  with  that 
of  a  geometrically  imperfect  spherical  shell: 


;u/2 

if- 033 


where  Ax  is  the  invariiint  spacing  of  stations  in  the  numerical  inte¬ 
gration. 

The  results  of  the  calculation  may  be  stated  as  follows  for  the 
critical  pressure 

where  -&z  is  the  local  maximum  pressure  at  which  buckling  occurs, 
^  is  the  buckling  pressure  of  a  complete  sphere  without  bending. 

For  a  perfect  spherical  shell,  one  finds  that 

In  terms  of  stresses,  one  can  write  that 


C & )  Perfect  shell 

*  ‘  '< ' 

O.-fSE  f-t\ 

=  /  3{i~VlV  Imperfect  shell 

This  means  that  the  predicted  buckling  stress  band  on  an  "imperfect 
shell"  analysis  is  approximately  80%  of  that  given  by  the  "perfect  shell". 
This  reduction  is  not  sufficient  to  bring  the  analytical  results  within 
the  range  of  available  experimental  results  for  elastic  hemispherical 
shells. 

However,  a  reduction  of  the  analytical  stresses  obtained  in  the 
present  viscoelastic  analysis  by  a  similar  magnitude  gives  a  more 
favorable  comparison  of  analytical  and  experimental  results  for  the 
hemispherical  shell. 
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Fig.  A.  3-1 

SHELL  GEOMETRY 
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EXPERIMENTAL  INVESTIGATION 
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II 


B.l  Experimental  Facility 


(a)  Heating  System 

The  functional  arrangement  of  the  experimental  facility  is  represented 
by  the  block  diagram  shown  in  Fig,  B-l. 

The  radiant  heating  oven  consists  of  a  polished  aluminum  reflector 
with  attached  heating  lamps.  The  heating  lamps  are  the  GE  T3/CL  quartz 
type.  Thirty- two  2000  T3/CL  lamps  form  the  heating  unit  for  Zone  I; 

Zones  II  and  III  each  consist  of  fifty-six  1000  T3/CL  lamps  (See 
Fig.  B-2  for  heating  zones). 

The  maximum  power  dissipation  of  the  entire  unit  is  478  KVA  which 
can  be  increased  by  replacing  the  lamps  of  Zones  II  and  III  with  2000 
T3/CL  lamps.  The  capacity  of  the  unit  corresponds  to  100  BTU/sq.  ft/sec. 
The  heat  input  to  the  model  is  less,  depending  upon  a  number  of  parameters 
including  distance  from  model,  model  configuration,  model  surface,  etc. 

The  efficiency  is  estimated  to  be  50-85%. 

The  oven  is  equipped  with  polished  aluminum  heat  shields  to  maintain 
zonal  temperature  control.  Two  sets  are  provided;  one  set  for  conical 
shells  and  the  other  for  hemispherical  shells. 

The  model  SPG6266W  three-phase  temperature  controller  provides 
zonal  control  of  temperature  in  the  oven  by  proportioning  the  a-c 
heater  power  through  ignitrons  controlled  by  a  closed-loop  thermocouple 
feedback  circuit.  Temperatures  may  be  automatically  maintained  to  within 
0.5%  by  employing  precision  potentiometric  circuitry  held  in  continuous 
calibration  by  a  standard  cell.  The  load  power  may  also  be  manually 
varied  from  zero  to  full  power  by  a  front  panel  dial.  The  temperature 
controller  is  used  to  program  temperature  as  a  function  of  time  by  switch¬ 
ing  to  the  remote  temperature  programmer.  The  temperature  controller 
can  proportion  a  continuous  load  power  of  390  KVA  for  the  three  phases 
with  water-cooled  ignitrons.  Each  load  phase  is  controlled  independent¬ 
ly,  thereby  providing  three  separate  zones. 

The  controller  is  a  complete  three-phase  ignitron  powered  control 
unit  combined  with  three  d-c  potentiometric  circuits  and  high  gain  a-c 
amplifiers.  Accurate  specimen  temperature  control  is  accomplished  by 
comparing  a  low  level  d-c  potentiometer  milli-voltage  with  the  thermal 
electromotive  force  from  the  specimen  thermocouple.  Any  positive  error 
between  these  voltages  is  chopped  and  a-c  amplified  to  a  relatively  high 
voltage  level.  The  amplified  a-c  output  voltage  is  converted  back  to 
d-c  in  a  diode  modulator  circuit  and  used  to  control  the  firing  angle 
of  the  ignitrons.  The  ignitrons,  in  turn,  proportion  the  power  to  the 
oven  to  restore  the  termocouple  to  the  set  point  temperature  and  to 
reduce  the  low  level  error  to  within  the  proportional  band  of  control. 
Compensation  for  ambient  temperature  chances  at  the  thermocouple  cold 
junction  is  automatic. 
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(b)  Temperature  Programmer 

The  Research,  Inc.,  Data-Trak  function  generator  (FGE5035)  provides 
an  output  proportional  to  the  ordinate  of  the  chosen  temperature -time 
function.  It  is  a  curve  follower  in  that  two  conducting  lines  bordering 
the  desired  curve  generate  a  high  electrostatic  potential  field  so  a 
null-seeking  probe  can  be  servo-driven  to  follow  the  zero  potential  line. 

The  temperature -time  graph  is  fastened  to  each  of  three  drums  such 
that  the  abscissa  (time)  is  the  rotational  position  of  the  drum.  The 
probe  is  mounted  between  two  parallel  slide  bars  beneath  the  drum.  As 
the  servo  positions  the  probe,  a  flexible  cable  from  the  probe  carriage 
rotates  a  precision  potentiometer  to  provide  an  output  voltage  propor¬ 
tional  to  the  translational  position  of  the  probe. 

Since  there  is  no  contact  between  the  conducting  lines  or  the  probe 
the  system  does  not  rely  upon  sensing  current  or  a  magnetic  flux. 

One  can  also  regulate  the  drum  cycle  from  10  to  300  seconds  pro¬ 
viding  a  flexible  time  scale.  In  addition,  by  placing  the  generator  in  . 
a  "hold"  condition  a  given  situation  can  be  maintained  for  as  long  as 
desired. 

(c)  Vacuum-Pressure  System 

The  vacuum  system  is  automatic  with  a  capability  of  permitting  a 
pressure  difference  across  the  shell  to  vary  within  tl/8  inch  of  mer¬ 
cury.  The  mercury  column  is  positioned  by  two  pins,  1/4  inch  apart, 
through  the  relay  and  the  vacuum  pump.  The  tank  acts  as  a  reservoir 
and  a  check  valve  while  the  relief  valve  allows  the  system  to  be  opened 
to  the  atmosphere  at  will. 

B.2  Instrumentation 


Modified  HT-600  strain  gages,  manufactured  by  the  High  Temperature 
Instruments  Corporation  are  used.  The  gage  is  modified  in  that  in 
addition  to  the  HT-600  alloy,  which  has  a  negative  coefficient  of  resis¬ 
tivity  with  temperature,  an  alloy  with  a  positive  coefficient  is  employed, 
thereby  reducing  the  apparent  strain.  The  gages  are  applied  in  a  mer¬ 
idional  and  parallel  direction  at  four  stations  along  a  meridian  using 
Allen  PBX  cement  (See  Fig  B-2).  Obtaining  a  good  bond  between  either 
the  pre-coat  and  the  aluminum  surface  or  the  second  coat  and  the  pre¬ 
coat  was  not  always  successful.  As  available  testing  time  became  short 
this  dictated  operation  with  less  than  full  strain  gage  instrumenta¬ 
tion.  Each  gage  used  a  three  wire  connection  to  the  bridge  thereby 
negating  effects  of  change  in  lead  wire  resistance.  The  signal  is  picked 
off  a  one  gage  bridge,  recorded,  and  then  corrected  for  temperature 
effects.  Since  the  apparent  strain  versus  temperature  curves  varied 
widely  for  different  gages,  each  gage  was  provided  with  individual  curves. 
Also,  the  nature  of  the  curves  changed  with  each  temperature  cycle  there¬ 
by  necessitating  strain-temperature  data  not  more  than  one  cycle  distant 
from  the  test  run. 
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To  calibrate  the  strain  gage  circuits,  a  decade  resistance  is 
shunted  around  the  active  gage  and  the  resistance  varied  while  the 
galvanometer  deflection  is  recorded.  Knowing  the  gage  factor,  one  can 
interpret  the  resistance  change  as  a  strain  and  the  deflection  then 
calibrated  in  terms  of  micro-inches/inch/inch  of  deflection. 

The  temperature  recording  circuit  is  simply  a  measurement  of  the 
current  induced  by  the  potential  difference  between  the  hot  and  cold 
junctions.  On  later  models  differential  thermocouples  (Fig  B-3)  provided 
data  on  the  temperature  difference  across  the  shell  wall.  One  junction 
is  a  spot  welded  to  the  Inner  surface  while  the  other,  after  being  led 
through  a  small  hole  in  the  wall  a  distance  from  the  measuring  area, 
was  cemented  opposite  to  the  previous  junction.  Since  one  junction  had 
to  be  electrically  insulated  from  the  other,  yet  in  good  thermal  contact 
with  the  outer  surface,  the  mechanical  bond  of  the  copper  oxide  cement 
used  often  failed, providing  only  qualitative  data  on  the  temperature 
differential.  An  aluminum  foil -backed  tape  shielded  the  outer  junction 
from  direct  radiation  while  also  sealing  the  hole  in  the  shell  wall. 

The  calibration  procedure  is  as  follows:  the  hot  and  cold  junctions 
are  brought  to  the  same  temperature;  a  voltage  is  induced  across  the 
thermocouple  resistance  and  the  galvanometer  by  a  precision  millivolt 
potentiometer;  with  the  aid  of  standard  potential  difference- temperature 
charts  for  the  particular  thermocouple,  the  circuit  is  calibrated  for 
°F  versus  deflection. 

All  the  instrumentation  wiring  in  the  model  was  supported  by  a  brass 
center  post  bolted  to  the  apex  of  the  model.  Brass  rod  branches  from 
the  center  post  supported  leads  to  the  individual  gages. 

A  porcelain  filter  disk  provided  a  means  of  passing  the  wiring  through 
a  pressure  differential  at  high  temperatures.  The  individual  leads 
passed  through  the  holes  in  the  disk  with  the  PBX  cement  acting  as  a 
potting  agent.  The  disk  was  then  bolted  to  the  steel  mounting  plate 
atop  a  quad  ring  providing  a  pressure  seal. 

B . 3  Test  Spec imens 

The  test  specimens  are  shown  in  Figs  B-4  and  B-5.  The  specimen  is 
attached  to  a  5/8  inch  thick  steel  plate.  A  silicone  rubber  "Quad"  ring 
is  inserted  around  the  specimen  support  to  provide  an  adequate  vacuum 
seal . 


The  specimens  are  coated  with  soot  to  increase  the  absorptivity  and 
to  provide  a  common  surface  for  all  tests. 

B.4  Material  Properties 

Although  the  model  manufacturer  and  Alcoa  assured  the  investigators 
that  soaking  the  models  for  a  period  of  thirty  minutes  at  650°F  anneals 
the  material  and  removes  any  residual  stresses  (produced  in  spinning),  a 
series  of  supplementary  tests  were  conducted  to  check  certain  important 
mechanical  properties  of  the  material. 

Models  were  tested  before  and  after  soaking  at  room  temperature 
under  the  action  of  a  concentrated  load.  It  was  felt  that  this  would 
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produce  the  desired  biaxial  stress  condition  without  complicating  the 
test  program. 

The  tests  were  limited  to  the  determination  of  the  modulus  of 
elasticity  and  of  Poisson's  Ratio.  Because  of  the  considerable  scatter 
of  test  data,  the  investigators  feel  that  these  data  are  valueless  to 
the  investigation  and  are  omitted. 


The  mechanical  properties  involved  in  the  analysis  are:  shear 
modulus  of  elasticity,  Poisson's  ratio,  coefficient  of  viscosity,  and 
the  coefficient  of  linear  thermal  expansion. 

(a)  Shear  Modulus  of  Elasticity 

The  following  approximate  values  for  shear  modulus  are  given: 
Temperature  (°F)  75  212  300  400  500  600  700 

G  x  10‘6  psi  3.8  3.7  3.6  3.4  3.0  2.6  1.5 


(b)  Poisson' s  Ratio 

At  elevated  temperatures,  Poisson's  ratio  is  a  multivaried  function 
of  elastic  strain,  plastic  strain,  creep  strain  (time  dependent),  and 
anisotropy  (varies  with  strain).  Poisson's  ratio,  therefore,  is  no 
longer  a  unique  property.  Tests  conducted  by  the  investigators  show  a 
range  of  values  from  0.33  to  0.45.  At  elevated  temperatures,  it  seems 
that  the  results  are  highly  sensitive  to  the  specimen  geometry  and  test 
procedure,  rather  than  a  strong  dependence  on  the  properties  of  the 
material . 

An  attempt  to  eliminate  such  variables  is  useless  without  an  exhaus¬ 
tive  number  of  well  controlled  tests. 

(c)  Coefficient  of  Viscosity 

No  attempt  was  made  to  establish  experimentally  the  coefficient  of 
viscosity  of  the  subject  material.  To  the  knowledge  of  the  writers, 
no  data  exist  of  the  variation  of  viscosity  with  temperature  below  the 
melting  point  for  aluminum  alloys.  Although  some  attempts  have  been 
made  to  establish  such  a  viscosity-temperature  curve  for  pure  aluminum, 
the  use  of  the  results  for  the  present  application  is  highly  questionable. 

It  is  agreed,  however,  that  the  viscosity-temperature  curve  for 
pure  metal  may  be  quite  different  from  that  of  alloys.  The  alloy  may 
show: 

1.  Different  rate  of  decrease  of  viscosity  with  increase  of  tem¬ 
perature. 
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2.  Different  rate  at  which  the  slope  of  the  viscosity- temperature 
curve  increases  as  the  liquid  state  is  approached. 

3.  A  greater  coefficient  of  viscosity. 

Based  on  tests  conducted  above  the  melting  points  (6),  the  results 
of  pure  aluminum  and  aluminum  alloys  indicate  the  following: 

1.  In  all  cases,  the  viscosity- temperature  curve  of  the  low  per 
cent  alloys  (9 aluminum)  is  similar  to  that  of  the  pure  metal. 

2.  For  higher  per  cent  alloys  (98%  or  less  aluminum),  the 
presence  of  the  alloy  had  a  characteristic  effect. 

Since  the  coefficient  of  viscosity  is  one  of  the  principal  parameters 
in  the  temperature-dependent  viscoelastic  analysis,  an  attempt  is  made 
to  estimate  the  viscosity-temperature  curve  for  the  material  of  the 
present  investigation.  An  expression  of  the  form 

>1  -  CG  exp(H/RT) 

1  2 

Where  C  is  a  constant,  G  is  modulus  of  rigidity  (dynes/cm  ),  H  is  heat 
of  activation  (calories/mole),  R  is  the  gas  constant,  and  T  is  the 
absolute  temperature  (°K),  is  used  to  represent  the  curve.  This  func¬ 
tional  relationship  is  a  typical  form  of  the  variation  of  viscosity  with 
temperature  for  some  polycrystalline  solids  and  liquids. 


In  the  present  investigation,  the  empirical  expression  was  taken  as: 
1  =  1.819  x  10”20  G  exp  (18000 /T) 

the  above  expression  is  based  on  established  values  by  Alcoa  (2)  at  the 
temperature  range  for  aluminum  alloys. 


The  graph  of  this  expression  is  shown  in  Fig.  2.1. 3-1  together  with 
the  results  obtained  by  Ke  (l6)  for  pure  aluminum. 

(d)  Coefficient  of  Thermal  Expansion 


No  attempt  was  made  to  investigate  the  coefficient  of  thermal  ex¬ 
pansion.  The  values  used  in  the  analysis  are  based  on  Reference  2  and 
communication  with  Alcoa. 


Fig.  2.1. 3-2  shows  the  variation  of  these  properties  with  temperature 
for  6061  aluminum  alloy. 
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TESTING  FACILITY  BLOCK  DIAGRAM 
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CONICAL  MODEL 


FIG.  B-2A 
INSTRUMENTATION 


110 


/TWO  'STRAIN  GAGE S 
/AND  ONE  THERMOCOUPLE 
'at  EACH  STATION 


Fig.  B-2B 

INSTRUMENTATION 


SHELL  WALL 
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Fig.  B-6  Instrumented  Conical  Shell  Model 


Fig.  B-7  Experimental  Facilities 


LEGEND  1.  Radiant  Heating  Oven,  2.  Power  Controller, 

3.  Temperature  Programmer,  4.  Vacuum  Controller,  5.  Zonal 
Heat  Shields,  6.  Oscillograph,  7.  Bridge  Balance  Unit, 

8.  Model. 
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TABLE  C-l 


MODEL 


TEADY  STATE  TEMPERATURE 


400°F 

400°F 

500°F 


(1)  400  (2)  300  (3)  200 


300°F 
400°F 
500°F 
400°F 
500°F  * 


(1)  500  (2)  400  (3)  300 


300°F 

400°F 

500°F 

400°F 

500°F 

(1)  400  (2)  300  (3)  200 
(1)  400  (2)  300  (3)  200 
(1)  500  (2)  400  (3)  300 
1)  500  (2)  400  (3)  300 
(1)  500  (2)  400  (3)  300 


1/16"  Hemispher 


(1)  400  (2)  300  (3)  20( 
(1)  500  (2)  400  (3)  30( 
(1)  500  (2)  400  (3)  30( 
.£11_5Q0_£2}_499_£3}_39( 


Repeated 


TEST  PROGRAM 


RATE  (°F/sec)  PRESSURE 


5  None 

10  None 

20  None 

(1)  10  (2)  5  (3)  2.5  None 


2 . 5  None 

5  None 

10  None 

10  22"  Hg 

20  22"  Hg 

(1)  50  (2)  40  (3)  30  22"  Hg 


2 , 5  None 

5  None 

10  None 

10  24"  Hg 

20  23"  Hg 

(1)  10  (2)  5  (3)  2.5  24"  Hg 

(1)  20  (2)  10  (3)  5  23"  Hg 

(1)  10  (2)  5  (3)  2.5  23"  Hg 

(1)  20  (2)  10  (3)  5  23"  Hg 

(1)  50  (2)  40  (3)  30  23"  Hg 


5  None 

5  21"  Hg 

10  None 

(1)  10  (2)  5  (3)  2.5  21"  Hg 

(1)  10  (2)  5  (3)  2.5  21"  Hg 

(1)  20  (2)  10  (3)  5  21"  Hg 

in_5Q_£2}_49_£3}_39J___18'_,_HS 


Temperatures  and  Stresses  in  l/32  Thick  Cone  for  a  Program  of  5  F/sec.  to  U00  F. 
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TABLE  C  -  3 

Temperatures  and  Stresses  in  l/32”  Thick  Cone  for  a  Program  of  10  Vsec.  to  U00  F 
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Temperatures  and  Stresses  in  l/32"  Thick  Cone  for  a  Program  of  20°  F/sec.  to  $00° 
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Temperatures  and  Stresses  in  l/l6  Thick  Cone  for  a  Program  of 
20  °F/sec.  to  500 °F  with  22"  Hg  Differential  Pressure 
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Temperatures  and  Stresses  in  l/l6  Thick  Cone  for  a  Program 
of  20  °F/sec.  to  fj00°F  with  22  Hg.  Differential  Pressure,  with  collapse 
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TABLE  C  -  20 

Temperatures  and  Stresses  in  l/8"  Thick  Cone  for  a  Zonal  Program  with  23"  Hg  Differential  Pressure 
Zone  1,  10  °F/sec.  to  500°F;  Zone- 2,  5  °F/sec.  to  UOO°F;  Zone  3,  2.5  F/sec.  to  300  F. 
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t  (  (  |HHHfvJcvjcvjcvjcvJcvjcvjHHHriHHHHH  ^  ,  ,  , 

1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

CVIQCOOvOOQOOOO-OOQOOOOOOOOOOOQQ 
vOvOnVLrVLTVCVJO  Win  CK4  H  O  vO  O  H,_c4  v&  co  H  win  H  O'.  cvvO  vO 
cvj  UV  vO  CO  O  _4  CO  H1TVO  CVI  HHco^JtM^vOvOlnJJtnrimm 
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ABSTRACT 


A  variational  equation  based  on  energy  principles  has  been  devel¬ 
oped  for  thin  hemispherical  shells.  The  behavior  of  the  shell  is 
assumed  to  be  characterized  by  a  viscoelastic  material  with  temperature 
dependent  properties. 

The  Ritz  method  is  applied  by  assuming  that  the  shell  displace¬ 
ments  are  functions  of  the  curvilinear  shell  coordinates  and  that  their 
unknown  coefficients  are  functions  of  time. 

A  system  of  six  ordinary  differential  equations  of  order  two  is 
developed.  These  equations  result  from  operations  on  the  varia¬ 

tional  equation. 
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D.l  Introduction 


A  number  of  investigations  have  been  directed  towards  determining 
the  thermoelastic  behavior  of  shell  structures.  In  recent  years  sev¬ 
eral  investigators^* 2, 3, 4, 5, 8, 9  have  used  variational  methods  to  treat 
problems  of  thermoelasticity  and  certain  linear  thermoviscoelastic 
problems . 


This  report  extends  some  of  these  results  toward  determining  the 
viscoelastic  behavior  of  shells  of  revolution  under  uniform  pressure 
and  thermal  gradients.  Although  certain  restrictions  are  made  concern¬ 
ing  the  coupling  of  heat  conduction  and  the  viscoelastic  phenomena,  the 
analysis  has  general  applicability. 


The  use  of  generalized  coordinates  and  variational  calculus  leads 
to  concepts  of  generalized  forces  of  the  Lagrangian  type  which  are 
applicable  to  mechanical  (static  and  dynamic)  and  thermal  (xsrith  material 
property  variations)  problems. 

First,  the  general  three -dimensional  thermoviscoelasticity  formu¬ 
lation  of  the  variational  equation  is  presented  in  curvilinear  coordi¬ 
nates.  The  energy  and  work  integrals,  including  the  material  properties, 
are  expanded  into  powers  of  the  normal  coordinate.  This  permits  the 
introduction  of  thermal  effects  on  the  material  properties. 


The  general  equations  are  specialized  to  thin  shells.  For  the  thin 
shell  analysis,  the  shell  thickness  is  allowed  to  approach  zero  and 
only  terms  of  shell  thickness  of  unit  power  are  retained.  An  example 
of  a  hemispherical  shell  is  given  to  illustrate  the  procedure  of  ana¬ 
lysis  . 


D . 2  Mathemat ic  al  Formul at ion 
D.2.1  Variational  Equation 


The  transient  response  of  a  shell  to  mechanical  and  thermal  loadings 
can  be  investigated  by  standard  energy  methods  without  the  complications 
connected  with  the  analytical  solution  of  the  partial  differential  equa¬ 
tions  of  equilibrium. 

If  the  shell  is  subjected  to  simultaneous  mechanical  and  thermal 
loading,  according  to  Hamilton's  principle  the  energy  formulation  is 
given  by  the  variational  equation 

(K+F-W)Jt  (D.2.1-1) 

Jio 
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where 


/  is  the  first  variation  and  integration  with  respect  to  time 

oerformed  between  the  limits  t  and  t, 

*•  o  1 

K  is  kinetic  energy 

F  is  free  energy  (potential  energy  in  isothermal  analysis) 

W  is  work  integral 

The  function  (K  +  F  -  W)  is  identical  to  the  Lagrangian  function. 
In  the  following  paragraphs,  the  proper  form  of  equation  D. 2.1-1  is 
developed  for  application  to  the  viscoelastic  analysis  of  shell  struc¬ 
ture  s . 


D, 2,1.1  Energy  and  Work  Integrals 

The  coupling  between  thermo elasticity  and  dynamics  is  expressed  by 
the  inclusion  of  the  acceleration  term  in  equilibrium  equation,  thus 
forming  the  usual  equations  of  motion 

C Ty  :  =  y*  it/  (D .  2 . 1  -  2) 

where  Ocj  and  Uc  are  the  unknown  field  variables.  This  leads  to  the 
Lagrangian  formulation,  expressed  in  terms  of  kinetic  energy 

f(^\/  (D. 2.1-3) 

where  f  is  the  density  of  the  material  and  represents  the  mass  of 

the  element  which  independent  of  time  ~t,  Hi  is  the  covarient  displacement 
vector,  and  (j/j  is  the  stress  tensor,  the  dot  represents  differentiation 
with  respect  to  time  t.  It  is  assumed  that  the  variation  of  the  density 
with  time  is  negligible. 


The  Helmholtz  free  energy  expression  to  be  used  as  the  potential 
energy  form  in  the  variational  formulation  is 


(D. 2.1-4) 


where  f  is  the  energy  density  per  unit  volume.  Equation  D. 2.1-4  can  be 
expressed  in  terms  of  the  stress  and  strain  tensors  as  follows 


F~  (D. 2.1-5) 

Now  let  us  introduce  the  stress-strain  relations.  The  material  is 
assumed  to  that  of  a  Kelvin  body  with  temperature  dependent  properties, 
in  which  case  the  following  function  for  the  stress  tensor  is  chosen. 
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(D. 2.1-6) 


Sty  tt  2  Cf&j  1 *  'j 
Ou=3*$t  +  3f£/f  j 

where  at  is  linear  thermal  expansion  coefficient,  1  is  coefficient  of 
viscosity,  G  is  shear  modulus,  K  is  bulk  modulus.  The  strain-dis¬ 
placement  relations  are 

Grifc+fy)  A"° 

Substitution  of  the  stress-strain  and  strain  displacement  relations 
into  the  free  energy  expression  yields 


(D. 2.1-7) 


To  include  the  effects  of  static  or  transient  external  loading, 
it  is  necessary  to  express  the  work  done  on  the  mass  of  material  by 
the  body  forces.  The  work  integral  may  be  written  as 


W-  f  fFc  (D. 2.1-8) 

where  F.  is  the  body  force  per  unit  m$ss  acting  on  the  material.  Or, 
in  terms  of  the  displacement  vector  u1  and/the  unit  pressure  intensity, 
we  find 


( P-uZAV  (D. 2.1-9) 

Jv  1 

It  is  assumed  that  the  only  external  force  acting  on  the  shell 
is  that  of  external  normal  pressure,  Pg.  Therefore  the  work  may  be 
written  as  a  two-dimensional  integral ,  where  now  P-,  is  a  function  of 
the  curvilinear  coordinates  y1  and  yz.  J 


W=  As  (D. 2. 1-10) 

Combining  equations  CD. 2. 1-4) ,  (D2.1-8)  and  (D2.1-10)  we  may  now 
write  equation  (D2.1-1)  as  follows 
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D.2.2  Geometrical  Relations  in  Curvilinear  Coordinates 


General  formulation  of  the  variational  equation  appropriate  to  shell 
theory  expressed  in  curvilinear  yi  coordinate  system  requires  the 
characterization  of  the  middle- surface  geometry.  Let  y*  represent  an 
orthogonal  curvilinear  coordinates  which  are  single-valued  functions 
of  the  Cartesian  coordinates  x1 

(D.  2.2-1) 

We  may  represent  any  point  on  the  mean-surface  by  yly2,  with  y3  as  the 
local  normal  as  shown  in  Figure  D.2.1. 

The  volume  element,  dV,  in  the  y1  coordinate  system  is  given  by 


(D.2.2-2) 

and  the  element  of  area,  dS  on  the  yi  surface  is 

4'42  (D.2.2-3) 

where  ylfyjl  is  the  square  roof  of  determinant  of  the  metric  tensor. 

In  equation  (2.2-2)^.  =  ^.  (yA,y  »y3)  while  in  (2.2-3)  ^  =  ^(y1,yZ). 


For  specific  examples  j/ has  taken  the  following  forms: 

(D.2.2-4) 
(D.2.2-5) 


(IU- 


A2  S/jJ  2 


Hemispherical  Shells: 

Conical  Shells: 

(vertical  height  is  unity)  where  Y  is  one-half  the  apex  angle. 
D.2.3  Deformation  of  a  Shell  and  Its  Middle  Surface 


=  -5/a'2  Y 


The  expressions  for  the  components  of  displacement  of  an  arbitrary 
point  of  the  shell  in  terms  of  displacement  of  the  corresponding  point 
of  the  middle  surface  are 


U 


i 


y\,  - 


«e'+u,  y 

//  2  » 

T  *1  / 


<c3= 


(D.2.3-1) 


•  .  ct  19 

where  Uq  and  ui  are  functions  of  yL  and  y4  only.  The  strain  displace¬ 
ment.  relations  in  curvilinear  coordinates  can  be  expressed  as 
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wnere 


f7{j  convarient  components  of  strain  tensor 

**  convarient  differential  coefficient  of  contravaxient  dis- 
placement  vector  f*  with  respect  to  the  curvilinear 
coordinates  y1 

metric  tensor 

Making  use  of  the  assumption  that  the  normals  to  middle  surface 
before  bending  remain  normal  after  bending,  we  get 


W*  *  2  t«) 


(D.2.3-3) 


After  performing  covarient  differentiation,  simplification,  and  sub¬ 
stituting  the  following  expression 


or  f  ' 

4^  -  H  j  H** 


(D. 2.3-4) 


we  obtain 


where 


/ 


7 


^bcT1^  '  £  "4 


—  —  2Z 


oC»4. 


(D.2.3-5) 


(D.2.3-6) 
no  summation  on  oi 


and 


L-t.  U 
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D.3  Variational  Equation  in  Terms  of  the  Displacement  of  the  Middle 
Surface  of  the  Shell. 

D.3.1  Expression  of  General  Displacements  in  Terms  of  the  Displacements 
of  the  Middle  Surface 


Consider  now  the  general  displacements.  Employing  the  results 
of  paragraph  D.2.3,  we  can  write  that 

■  v  '  '  u/  '  '■  7  i  '  ('  j 

-  M}]+ 

*  IiT\e  (l;  iijj + 

+  E“('lyJ+2fE“  ISfU  «/+  (C  “.‘j  J  (d.  3.1-1) 

o-3-1-2) 

From  equation  (J). 2. 3-5)  we  have 

K-‘C-f({Fij£*G'C) 

«3-«.3 

Differentiating  and  substituting  above  result  we  can  express 
equation  (D3.1-2)  as  follows: 

vv  *  (% -itf*  (cT+ 

(D. 3.1-3) 


where 
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C=*  ft®  I  (D.  3.1-4) 

//  =  j/£”  |  (D.3.1-5) 

D.3.2  Variational  Equation  for  a  General  Shell  of  Revolution 


By  making  use  of  equations  (D. 2. 1-11),  (D.2.2-2),  (D.3.1-1)  and 
(D. 3.1-3),  we  can  write  the  variational  equation  for  a  general  shell 
of  revolution  (shell  of  general  geometry)  in  the  following  form: 


(D.  3.  2-1) 
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D.3.3  Variational  Equation  for  a  Hemispherical  Shell 


Now  we  expand  equation  (D.3.2-1)  for  application  to  hemispherical 
shells  in  terms  of  middle  surface  deformations. 

From  Equations  (D. 2.3-1)  and  (D.2.3-5)  we  have  that 


XO 


-  2  L 


AA 


9*«  I 

\f,o 


(D.3.3-1) 


£  =  /  s  -SluS 

L-"  "Ifio 
/ 


£  •  4-  =  -V 


E  =  4.,./*  —  -  JL  Z 


®  -a 


I  (D.3.3-2) 


From  the  above  two  sets  of  equations  we  can  write  that 


(D. 3.3-3) 


(D.3.3-4) 


The  functions  A,  B,  C,  D,  F  and  H  take  the  following  forms: 
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AT  so  for  a  hemispherical  shell 


After  putting  these  values  in  (D.3.2-1),  simplifying  the  expres¬ 
sion  and  integrating  with  respect  to  y3  through  the  thickness  of  the 
shell,  we  get  the  following  variational  equation  for  a  hemispherical 
shell. 


Variational  equations  for  the  hemispherical  shell 
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where 


Ts~  U)  is  surface  temperature 


D.4  Integration  of  the  Variational  Equation  for  Thin  Hemispherical 
Shell 


D.4.1  Boundary  Conditions 


Consider  a  hemispherical  shell  which  is  clamped  at  the  edge,  0  = 
(See  Fig.  D.4.1).  At  0  ="%  -and  04&421T,  the  deformations  u,  v,  w 
and  the  rotation  ijr  become 

H-lTn  ou  =  o 

+  J£ 

In  case  of  a  hemispherical  shell  we  have  that 


«=  4'=  «,'+  -%-( V-£T/  T&) 


U>  -  LL  =  a,t  *  0  f  c<2  /i,  £ 

Az"  ^  ^  OK.  /?.,=  ''*• 

Considering  now  the  boundary  conditions,  we  can  write  that 
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which  can  be  written  as 


(E. 4.1-1) 

(D.4.1-2) 

(D.4.1-3) 

(D.4.1-4) 

From  equations  (D.4.1-2)  and  (D.4.1-4) 


Hence 


OR 


kc  =  « 


k*  -  o 


Therefore 
Similarly 

And  from  equations  (D.4.1-3)  k0  **° 

Thus  we  have  the  following  boundary  conditions: 

ko  -m  kp  „  u2  =  O 


Tjh- 


—  o 


f  at  AuD 


Also  the  functions  u^,Un,Uo  must  be  periodic  with  respect  to  9 
with  a  period  2  JT. 

D.4.2  Expression  of  Displacements  in  Terms  of  Functions  0  and  9 


It  is  necessary  to  represent  u.0  f  uj  ;  ua  in  terms  of  functions 
of  0  and  9  whose  coefficients  are  unknown  functions  of  t.  These  func¬ 
tions  must  satisfy  the  boundary  conditions  and  should  be  periodic  with 
respect  to  6  with  a  period  of  2 1C .  These  conditions  are  fulfilled  by 
the  following  functions: 


61 


'■ia  -u  Coi  ly?  Cos  $  S/*!1  & 

Cos  +  $2.  £aij  Si & 

<4  -  ( ?-¥){j h  coih£  “*/*"'*) 

These  functions  or  their  derivatives  are  substituted  in  the  var¬ 
iational  equation  for  a  thin  hemispherical  shell  and  integration  is 
performed  with  respect  to  @  and  G  over  the  given  boundaries. 


D.4.3  Integration  of  the  Variational  Equation  for  Thin  Hemispherical 
Shell 


For  the  hemispherical  shell  the  integration  has  been  carried  out 
over  the  limits  of  S  &  and.  .  But  as  we  must  take  the  limits 

for  as  e.  go  around  the  circle,  a  factor  of  4  will  be  common 

in  all  integrations.  Since  it  is  a  common  factor,  it  will  not  affect 
our  results.  Some  of  the  integrations  have  been  carried  out  numer¬ 
ically.  A  singularity  occurs  at  the  apex  and  in  one  of  the  integrations 
an  approximate  value  lias  been  taken  numerically  a  surface  making  an  angle 
of  <j>  =  0°30'  with  the  axis.  Thus  the  results  of  this  analysis  cannot  be 
applied  to  the  apex  and  a  separate  analysis  will  be  needed  fromthis  region 
of  the  shell. 

After  having  carried  out  the  integration  over  the  middle  surface 
of  the  shell,  we  reduce  the  variational  equation  to  functions  of  time 
only.  This  equation  is  shown  as  follows:  (See  next  page). 
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Variational  Equation  for  a  Thin  Hemispherical  Shell 

+  4 *6  j«- 2 Ufa + ia*?fa  /. *f7/Stfc *  *  3pfil +0.U/  fa  loftf  %  /  *4/ fa  +  6.  3 p  fa 

t  t  <fk I a t#fi+  *. r px+e. tjtfil  + 

t-  'itf-firs  *  ^ifas*  Tl*(*s 2 %£+£')+*&*  frsKfsrjsfify 

t  o.  if*  (fa  zt'x  ifxy  (o.vif,  v  «  mr,& + *•  *?/£)  *  *tft(r,  fa,  z) + 

+  *-w  (/&+*/:& &*;)+■  *2*2fefaA)f  +  a-i6/fas '■  °*ip,  fa0-  u<(p,fa 

+ (S,  pj + *.  ip  fas*  *‘i  fas  in 7  fas0**  (fas  tty*-  >v  KK+*  2??fJs 

+ °^r,t  *  °^(r,far,fz)  *  *-  “jnh .}+?*  h  **M /  aMz(^fafa-) + 

*  "?tfas‘t}7f> i  - “tyfak +ffl+  * ■#*{/*£+&*)+• *-rfKfrfi+Kf.Hi)+ 

+  sfs^yaH2(rfaiKy0sc(rX^ 

rfa)+  +• 

+  0-ty(rj+tth*-nt(r,rsfc)+ °-w  (r,rsf,A)+  o-nrfcrsM )  + 

r  *-X*f*A  -  '■  w/?*£+£rz)J  + 

+  /.  rp  rJ-  f^sctTs  i  r7/iyV'  ySTst  +  fk(f.  t&fi  +*  f*jijJ  Jt  =  o 


(D.4.3-1) 


where  fa  ft  (b-f)[/-  j(gfj 


Tg  Is  surface  temperature  and  function  of  t  only. 
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D.5.  Development  of  Different tial  Equations  for  Thin  Hemispherical  Shell 


D.5.1  Functional  of  the  Variational  Equation 


Tlie  variational  equation  can  be  represented  in  the  following  form: 


(D.5.1 =1 


The  integrand  F  is  known  as  functional.  It  consists  of  indepen¬ 
dent  variable  t  and  other  known  and  unknown  functions  of  t  e.g. 

n  f  etc.  Expression  for  the  functional  F  is  shown  as 
follows:0  ‘  ' 


Functional  F  of  the  Variational  Equation  is 


Fm  05i4fasj4irfeojfir*+oA7tj>  v ofytfr + J  + 

-t  fMk  ii/ fa  hfitf+iMjfip+i’-tyft-f  * 

+e.44fyso.n7A-f 

-wr.z-  ‘■xiffih*  mrstf)*  **■(&■/&)}+ 

-#■  ^Yaa+M)*  swf  i+'-vfKit+tj:)* 

+  »■  + *  Wforsfa)  r/i  +r,$+  i(z/2£  ~fiw(U+/’4)J  + 

+  !.q  1^4*7?+  +Fn(/-r&f/+0.f2£flj 


where 


are  functions  of  time 
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D.5.2  Euler  Lagrange  Equations  for  the  Thin  Hemispherical  Shell 


We  assume  that  value  of  functional  F  vanishes  at  both  the  time 
t  -  to  and  t  =  tj.  It  is  on  this  assumption  that  we  can  make  use  of 
Euler -Lagrange  equations  for  expression  of  the  results  in  form  of  dif¬ 
ferential  equations.  The  Euler -Lagrange  equations  for  the  functional  F 
with  unknowns  6.  f  and  their  time  derivatives  can  be  expressed 

in  the  following  form. 


(D.5.2-1) 

(D.5.2-2) 

(D.5.2-3) 

(D.5.2-4) 

(D.5.2-5) 

(D.5.2-6) 


The  Euler-Lagrange  Equations  are 

+hs'?+°nr)lA+ Unfa  -  ° 

4  for*  fefi+ij  6fYt~  3.3$^^-/.  UCQfe  -/■  (2  *  V  fyK* 

+(t+.os)nf,+  U-jit+OMfa  *  * 

-f  {2Snr+ojsjja^-h  U-Jif+  '-Zfjfi  -  0 

f£hf+0Sf/ihj>- 

_  l.fffii  -t (o.yi/f'-Hf.p'Jyht  +  A 4irL+0.ftJj / fa  + 

-f- /f3j  b^+O.  fZj  (D. 5. 2-11) 

orjC^riL'^  n+4-rtffif  -feXW-f  *  7^44  ~ 

~3So6ft>rr  2.j+4ttl+{o.f?rL+ &fi)f+(qo<?+ 

-3.(>fy  +  /SJjti' t/ljjkr2-3-Z/pn*o{ D.  5. 2-12) 


(D.5.2-7) 

(D.5.2-8) 

(D.5.2-9) 


.'r.  ^  **  *1  n\ 

■  J.  4-l.V/y 
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D.5.3  Discussion 


The  general  variational  equation  for  the  investigation  of  the  visco¬ 
elastic  behavior  of  shells  of  revolution  subjected  to  uniform  pressure 
and  thermal  gradients  are  developed.  The  development  is  based  on  the 
assumption  that  the  material  behaves  as  a  Kelvin  body  with  temperature 
dependent  properties.  The  equation  is  applied  to  a  shell  of  general 
geometry. 

By  making  use  of  the  geometrical  relations  for  a  hemispherical 
shell,  the  general  equation  is  specialized  for  application  to  thin  hemi¬ 
spherical  shells.  This  is  accomplished  by  neglecting  all  higher  order 
terms.  The  respective  geometrical  relations  for  conical  shells  are 
presented  for  a  similar  application. 

The  boundary  conditions  applied  in  the  investigation  are  based  on 
a  clamped  shell.  The  Ritz  method  is  applied  by  assuming  that  the  dis¬ 
placement  functions  which  satisfy  the  boundary  conditions  and  their 
respective  coefficients  are  dependent  on  time.  Other  boundary  conditions 
can  be  incorporated  by  an  appropriate  choice  of  the  displacement  func¬ 
tions. 

Integration  of  the  variational  equation  (now  a  function  of  the  dis¬ 
placements  and  their  respective  functions)  over  the  limits  governed  by 
the  shell  geometry,  results  in  an  equation  of  unknown  functions  of  t. 

Prom  the  variational  equation  a  system  of  Euler  Lagrange  equation 
is  developed.  In  the  specific  problem  we  have  a  system  of  six  ordinary 
differential  equations  of  order  two.  The  coefficients  of  the  unknown 
functions  are  various  functions  of  time.  The  equations  are  assumed  to 
be  continuous  and  that  their  values  vanish  at  t  =  tQ. 

A  solution  of  this  system  of  equations,  which  determines  the  un¬ 
known  functions,  must  be  made  numerically  by  digital  computer.  The 
sol  ution  is  based  on  the  time  interval  <<t^  tp 
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